!7flb-fil24  189 


UNCLnSSIFIEO 


FACTORS  INFLUENCING  THERHONECHANICAL  FAILURE  OF  FACE 
SEALS  II<U}  THAVER  SCHOOL  OF  ENGINEERING  HANOVER  N  H 
F  E  KENNEDV  ET  AL.  JAN  83  N8aei4-81-K-0e9a 

.  F/G  11/1 


1/1 


NL 


m 

END 

V-  *• 


FACTORS  INFLUENCING 

THERMOMECHANICAL  FAILURE  OF  FACE  SEALS  II 


Francis  E.  Kenne4yt  Jr. 
*#i<>4Assottate  Professor  of  Engineering 

r  ■  .... 


'.■■*  -  J:"  ‘'Vr,' 


‘‘iSfa^sto  l^pilawrh  ^^:V, 


. ; -Vv  ,.(»  *•^2 


Interim  Report  #2 
.submitted  to 


Office  of  Naval  Research 
Contract  No.  N00014-81-K-0090 
Period  Covered:  December  1,  1981  to  November  30,  1982 


DTIC 

EUECTE! 
FEB  8  19^ 


January  1983 


~DrSTRIiiUTl6N  STATEMENT  A 

iVppioved  !o«  public  trieoM( 
Difltflbutiott  Unlimited 


015 


Unclassified 

SKCUftiTV  classification  of  this  P«GE  fW7i«n  Pl«  Entmnd) 

I  REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NOMBER 


2.  COVT  ACCESSION  NO. 

/4/P^//  91 


4.  title  Cand  SubtUU) 

FACTORS  INFLUENCING  THERMOMECHANICAL 
FAILURE  OF  FACE  SEALS  II 


READ  INSTRUCTIONS 

_ BEFORE  COMPLETING  FORM 

3.  RECIPIENT'S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Interim  Report  #2 

Dec  K  1981  -  Nov.  30.  1982 

t.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORfa) 

Francis  E.  Kennedy,  Jr. 

John  N.  Grim 
Roger  P.  Glovsky 

s.  performing  organization  name  and  address 

Thayer  School  of  Engineering 
Dartmouth  College 
Hanover,  New  Hampshire  03755 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Engineering  Science  Directorate 
Office  of  Naval  Research 
Arlington,  Virginia  22217 

14.  MONITORING  AGENCY  NAME  A  ADDRESSfK  dIMarwM  Inm  CcntrolUng  Olltcm) 


A.  CONTRACT  OR  GRANT  NUMBERr*; 


N00014-81-K-0090 


to.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 


12.  REPORT  DATE 

January  1983 _ 

13.  NUMBER  OF  PAGES 

72 

IS.  SECURITY  CLASS,  (ol  tMm  nport) 

Unclassified 


I  It.  OISTRIBUTION  statement  (ol  IMt  Rtport) 


ISa.  DECLASSIFICATION/OOWNGRAOING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited. 


I  n.  OISTRIBUTION  STATEMENT  (ot  th0  mb9trmct  In  BfoeJk  20,  1/  dliUrmt  from  R^pptf) 


ELECTEIf 

FEB  8  1983  ’ 


IIS.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Conflnu*  on  rororoo  oldo  l/nocoo«orr  ond  Idonlliy  by  bloelt  numSor) 

Mechanical  Seals,  Face  Seals,  Seals,  Therrooelastic  Instability, 

Failure,  Surface  Temperature. 

.  4. _ _ _ 

Analytical  and  experimental  studies  of  factors  affecting  thermocracking 
and  other  modes  of  thermomechanical  failure  of  mechanical  face  seals  have 
continued  .for  a  second  year _ 

In  the  experimental  phase  'a  unique  contact  probe/ thermocouple  has  been 
built  for  studying  characteristics  of  small  patches  of  solid-to-solid  contact 
at  the  sealing  interface  of  mechanical  face  seals  during  seal  operation.  The 
probe  was  used  to  determine  the  nimber,  size,  location,  and  surface  y  (pto) 


DD  I  JAN ^3  1473  edition  OF  I  NOV  AS  IS  OBSOLETE 


S/H  0102-LP-OI4.A601 


Unclassified  J 

security  CLASSIFICATION  OF  THIS  PAGE  rWTiM  DM*  EnIWddJ 


Unclassified _ _ 

tgCU/<ITV  classification  of  this  page  flWiwi  Data  gnffrfj _ '  _ _ _ 

20.|Xcont1nued) 

temperature  of  contact  patches  over  a  range  of  seal  velocities.  Both  dry 
and  liquid-lubricated  conditions  were  studied,  and  three  different  mating 
ring  materials  were  tested.  It  was  found  that  small,  discrete  coi>’:act 
patches  are  indeed  present  on  the  contact  interface,  that  the  patches  tend 
to  be  stationary  with  respect  to  the  metallic  mating  ring,  whether  that  ring 
is  stationary  or  rotating,  and  that  the  size  and  temperature  of  the  contact 
spots  is  influenced  by  operating  velocity  and  by  the  thermal,  elastic,  and 
wear  properties  of  the  seal  materials.  . 

The  analytical  phase  of  the  study  was  focused  on  the  prediction  of 
surface  temperatures  at  the  patches  of  contact  on  the  seal  interface. 
finite  element  surface  temperature  analysis  program  (THERMAP)  was  rew^tten, 
documented,  and  proven  out  on  a  number  of  test  problems.  The  program  was 
used  to  study  the  influence  of  various  operating  and  material  parameters  on 
the  magnitude  of  the  surface  temperatures  in  the  contact  patches  of  face 
seals.  -Jt  was  found  that  the  most  promising  approach  to  decreasing  the 
surface  temperatures  would  be  to  increase  the  thermal  conductivity  of  the 
mating  ring  and  primary  ring  materials.  A  decrease  in  seal  velocity  (by  a 
decrease  in  diameter)  and  improved  cooling  of  the  seal  faces  would  also 
serve  to  decrease  surface  temperatures.  It  was  also  found  that  numerical 
oscillations  can  develop  in  the  solution  at  high  Peclet  numbers  (high 
velocity,  low  diffusivity,  large  finite  elements).  Several  schemes  were 
studied  for  the  elimination  of  the  oscillations  and  it  was  found  that  a  | 
streamline  upwind  procedure  is  quite  effective  in  most  cases.  j 


Unclassified 


SKCURITY  CLAfSIFICATION  OF  TMlI  RAOgfWi**!  Knfnd) 


FACTORS  INFLUENCING 

THERMOMECHANICAL  FAILURE  OF  FACE  SEALS  II 


Interim  Report  #2 
submitted  to 

Office  of  Naval  Research 
Contract  No.  N00014-81-K-0090 
Period  Covered:  December  1*  1981  to  November  30,  1982 


by 

Francis  E.  Kennedy,  Jr. 
Associate  Professor  of  Engineering 


and 


John  N.  Grim  and  Roger  P.  Glovsky 
Graduate  Research  Assistants 


Thayer  Scnool  of  Engineering 
Dartmouth  College 
Hanover,  New  Hampshire  03755 


January  1983 


Reproduction  In  whole  or  In  part  Is  permitted  for 
any  purpose  by  the  United  States  Government 


FOREWORD 


Work  at  the  Thayer  School  of  Engineering  at  Dartmouth  College  on  this 
project  has  been  supported  by  Office  of  Naval  Research  Contract  No.  N00014'81-K 
Mr.  M.K.  EHingsworth  of  ONR  is  the  contract  monitor. 

The  authors  gratefully  acknowledge  the  cooperation  of  personnel  from  the 
David  Taylor  Naval  Ship  R&D  Center,  especially  in  the  early  stages  of  the 
project.  In  particular  Mr.  Sidney  A.  Karpe  of  DTNSRDC  assisted  in  microscopic 
observations  and  contributed  to  many  fruitful  discussions  about  microscopic 
aspects  of  thermocracking. 

Mr.  Victor  A.  Surprenant  of  Dartmouth  College  assisted  in  metallography 

1 

and  in  interpreting  the' microscopic  observations. 

Face  seals  for  the  experimental  phase  of  the  project  were  contributed  by 
EG  &  G  Sealol,  Inc.,  H.F.  Greiner  Vice  President  and  Chief  Engineer. 

The  principal  investigator  (FEK)  gratefully  acknowledges  the  use  of  the 
computer  facilities  at  the  Laboratoire  de  Mgcanique  des  Contacts,  INSA,  Lyon, 
France  for  a  portion  of  the  computational  work  on  this  project. 


Accession  For 

NTIS  GRAScI 
DTIC  TAB 
Unannounced  D 

Justificotlon - 


-0090. 


TABLE  OF  CONTENTS 


INTRODUCTION  . 

SUMMARY  OF  PREVIOUS  RESULTS  ON  CONTRACT  N00014-81-K-0090  - 

EXPERIMENTAL  STUDY  OF  CONTACT  PATCHES  IN  FACE  SEALS  . 

EXPERIMENTAL  APPARATUS  . 

TEST  PROCEDURES  AND  MATERIALS  . . — - 

RESULTS  AND  DISCUSSION  OF  EXPERIMENTAL  STUDY  . 

Rotating  mating  ring-dry  - 

Rotating  carbon  primary  ring-dry  - 

Rotating  metallic  ring  with  sealed  fluid  - 

CONCLUSIONS  FROM  EXPERIMENTAL  STUDY  . 

ANALYTICAL  STUDY  OF  TEMPERATURES  NEAR  CONTACT  PATCHES  - 

FINITE  ELEMENT  THERMAL  ANALYSIS  PROGRAM  (THERMAP)  - 

FACTORS  INFLUENCING  SURFACE  TEMPERATURES  IN  FACE  SEALS 

EFFECT  OF  MOVING  SURFACE  VELOCITY  ON  THERMAL  ANALYSIS 

Evidence  of  numerical  Instabilities  at  high  velocity  - 

Techniques  for  handling  convective-diffusion  equations 

Application  of  upwinding  techniques  to  THERMAP  -  -  -  - 

CONCLUSIONS  FROM  ANALYTICAL  STUDY  . . 

REFERENCES  . - . - . 

APPENDIX  A  THEORETICAL  DEVELOPMENT  OF  EQUATIONS  - - 

APPENDIX  B  THERMAP  FLOWCHART  . . 

APPENDIX  C  USER  INSTRUCTIONS  FOR  THERMAP  . . 

APPENDIX  D  LIST  OF  PUBLICATIONS  AND  PRESENTATIONS  RESULTING 
FROM  THIS  SEARCH 

APPENDIX  E  DISTRIBUTION  LIST  . * 


LIST  OF  FIGURES 


Figure  Page 

1.  Schematic  diagram  of  contact  probe  (above)  7 

Idealized  output  signal  (below). 

2.  Photograph  of  stationary  carbon  primary  seal  ring  containing  contact  10 

probe,  in  its  holder.  Electrical  contact  brush,  which  slides  against 
mating  ring,  is  also  included. 

3.  Schematic  diagram  of  entire  test  system.  12 

4.  Comparison  between  contact  probe  output  (Fig.  4a)  and  appearance  of  16 

mating  ring  surface  (Fig.  4b).  Test  A,  440C  stainless  steel  mating 

ring,  188.5  rad/sec,  dry. 

5.  Comparison  between  contact  probe  output  (Fig.  5a)  and  appearance  of  16 

mating  ring  surface  (Fig.  5b).  Test  B,  440C  stainless  steel  mating 

vring,  188.5  rad/sec,  dry. 

6.  Comparison  between  contact  probe  output  (above)  and  thermocouple  18 

output  (below). 

Test  E,  440C  stainless  steel  mating  ring,  125.7  rad/sec,  dry. 

7.  Sequence  of  output  traces  from  contact  probe,  showing  growth  of  contact  21 
patches  to  steady  state  size. 

8.  Contact  probe  output  (steady  state).  22 

Beryllium  copper  mating  ring,  188.5  rad/sec.  Test  H. 

9.  Contact  probe  output  (steady  state).  22 

52100  bearing  steel  mating  ring,  188.5  rad/sec.  Test  0. 

10.  Sketch  of  primary  seal  ring,  showing  section  studied  in  thermal  and  30 

stress  analyses. 

11.  Finite  element  mesh  for  mechanical  face  seal.  31 

12.  Isotherms  in  mating  and  primary  seal  rings  for  case  of  Figure  11.  34 

13.  Study  of  velocity  effect  in  frictionally  heat  generating  contact  problem.  40 

14.  Finite  element  test  models  for  study  of  velocity  effect.  45 


LIST  OF  TABLES 


Table  No. 

1.  Properties  of  mating  ring  materials  used  in  this  study. 

2.  Characteristics  of  steady-state  contact  patches  for  various 
mating  ring  materials. 

3.  Parametric  study  for  mechanical  face  seal  analysis. 

4.  Results  of  thermal  analyses  of  moving  slab  subject  to 
distributed  heat  source. 

5.  Results  of  one  dimensional  thermal  analysis  with  moving  heat 
source. 

6.  Moving  slab  with  distributed  heat  source  -  results  of  analysi 
with  upwinding. 


1 


INTRODUCTION  AND  BACKGROUND 
Introduction 

Mechanical  face  seals  are  designed  to  prevent  fluid  leakage  around  rotat¬ 
ing  shafts  at  points  where  the  shaft  passes  through  a  stationary  housing  or 
pressure  vessel.  They  have  proven  to  be  the  most  effective  type  of  seal  for 
severe  applications  involving  high  fluid  pressures,  high  shaft  speeds,  or  a 
combination  of  these  conditions.  Unfortunately,  as  the  severity  of  the  appli¬ 
cation  increases,  the  frequency  of  seal  failure,  as  indicated  by  excessive 
leakage  and/or  excessive  seal  friction,  also  increases.  Such  failures  should 
be  avoided,  if  possible,  especially  in  marine  applications  where  the  seal  pre¬ 
vents  Intrusion  of  sea  water  along  a  rotating  propeller  shaft. 

Several  of  the  most  common  modes  of  face  seal  failure  are  caused  by  a 
combination  of  high  temperature  and  high  stress  in  the  seal  components.  These 
thermomechanical  failure  modes,  which  include  thermocracking,  or  heat  checking, 
and  warping  of  the  mating  ring,  are  often  responsible  for  excessive  seal  leak¬ 
age,  excessive  wear  of  the  primary  ring,  and,  occasionally,  fracture  of  the 
mating  ring.  Attempts  to  eliminate  thermocracking  and  related  thermomechanical 
problems  in  face  seals  have  not  been  completely  successful,  primarily  because 
our  understanding  of  these  mechanisms  is  very  limited.  The  goal  of  this 
research  has  been  to  gain  a  better  understanding  of  factors  affecting  thermo- 
cracking  and  other  related  failure  mechanisms  in  mechanical  face  seals. 

It  has  been  known  for  a  number  of  years  that  thermocracking  in  face  seals 
is  caused  by  thermal  stress  resulting  from  frictional  heating.  Several  studies 
of  the  problem  [1,2]  resulted  in  the  proposal  of  material  parameters  by  which 
a  material's  resistance  to  thermocracking  can  be  predicted.  Those  studies  did 
not  really  establish  a  mechanism  for  thermocracking,  however,  and  they  assumed 
that  temperature  and  stress  distributions  on  the  sliding  interface  are  both 
uniform  -  an  assumption  which  has  since  proven  to  be  false. 
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In  order  to  achieve  good  performance  from  a  mechanical  f^e  seal,  the 
sealing  surfaces  of  both  primary  seal  ring  and  mating  ring  are  lapped  to  a  high 
degree  of  flatness  during  production.  Despite  this,  a  very  slight  amount  of 
initial  waviness  remains  and  this  waviness  has  been  observed  to  grow  during  seal 
operation  [3].  One  consequence  of  this  waviness  is  that  frictional  surface 
tractions  are  not  uniformly  distributed  over  the  sealing  Interface  and  this 
results  in  circumferential  variations  in  frictional  heating  of  the  seal  rings. 

It  was  shown  by  Burton  and  his  co-workers  [4]  that  under  certain  conditions  this 
nonuniform  heating  may  lead  to  even  more  concentration  of  friction  and  frictional 
heating,  and  to  the  formation  of  hot,  highly  stressed  contact  patches,  or  thermal 
asperities.  Later  analytical  work  determined  the  influence  of  operating  vari¬ 
ables,  such  as  velocity,  and  of  material  parameters,  such  as  thermal  conductivity, 
modulus  of  elasticity,  coefficient  of  thermal  expansion,  wear  resistance  and 
friction  coefficient,  on  this  unstable  process  of  transition  from  nominally  flat 
to  highly  deformed  surfaces  [5-7J. 

These  analytical  studies  showed  conclusively  that  this  unstable  condition, 
called  thermoelastic  Instability,  could  occur  In  face  seals.  It  was  hypothesized 
that  the  stresses  and  temperatures  near  the  resulting  thermal  asperities  could 
be  responsible  for  thermocracking  and  other  modes  of  seal  failure.  There  was 
some  evidence,  based  on  observation  of  disassembled  face  seals  after  failure, 
that  this  was  In  fact  what  had  occurred  [8].  The  hypothesis  was  not  proven 
conclusively  because  several  things  were  missing:  first,  experimental  evidence 
that  small,  hot  contact  patches  resulting  from  thermoelastic  instability  actually 
occur  In  operating  face  seals  and.  If  so.  Information  about  factors  influencing 
the  size,  location,  and  number  of  these  patches  and  contact  conditions  within 
them;  and  second,  a  study  showing  how  conditions  near  the  contact  patches  can 
lead  to  the  onset  of  t’  -"ncra  ng.  This  research  project  has  been  attempting, 
with  some  success,  to  provide  this  Information  and  to  establish  a  firm  basis 
for  reconmendatlons  for  the  design  of  face  seals  which  would  not  be  subject  to 
a  high  Incidence  of  thermomechanlcal  failure. 


SUMMARY  OF  PREVIOUS  RESULTS  ON  CONTRACT  N00014-81-K-0090 

Since  December  1,  1980,  work  has  been  in  progress  at  Thayer  School  of 
Engineering  on  ONR  Contract  N00014-81-K-0090.  During  the  first  year,  the  research 
was  directed  at  learning  the  root  causes  of  thermocracking,  whether  these  causes 
are  related  to  non- uniformity  of  contact  at  the  seal  interface,  and  what  material 
parameters  have  the  most  influence  on  thermocracking  mechanisms.  The  research 
included  both  analytical  and  experimental  activities  aimed  at  accomplishing  these 
goals.  Most  of  the  results  achieved  during  the  first  year  of  the  research  program 
have  been  published  elsewhere  [9,10].  Those  results  will  be  summarized  here. 

The  first  phase  of  the  work,  accomplished  with  the  aid  of  personnel  at 
the  David  Taylor  Naval  Ship  RiD  Center  in  Annapolis,  involved  the  examination  of 
marine  seal  components  which  have  exhibited  thermocracking.  Macroscopic  and 
microscopic  examination  of  the  metallic  seal  rings  showed  that  most  thermocracks 
are  radial  in  nature  and  initiate  at  carbide  particles  on  the  seal  surface  cf  the 
cast  cobalt-based  alloy  rings  [9].  The  cracks,  which  are  distributed  relatively 
uniformly  in  the  circumferential  direction  on  the  seal  face,  appear  to  initiate 
and  propagate  in  a  brittle  manner  in  a  tensile  mode  (Mode  I).  There  was  evidence 
on  the  seal  surfaces  that  localized  contact  conditions  had  occurred  in  thermo- 
cracked  regions,  with  both  deformation  and  temperature  in  the  contact  region 
being  rather  severe. 

The  next  phase  of  this  research  was  an  analytical  study  which  aimed  to 
determine  how  localized  contact  patches  at  the  seal  interface  could  lead  to  the 
onset  of  thermocracking.  The  finite  element  method  was  chosen  and  a  finite 
element  thermal  analysis  program  was  developed  for  use  in  calculating  surface 
temperature  distributions  in  sliding  seals  [11].  This  program,  THERMAP,  was 
used  along  with  a  finite  element  stress  analysis  program,  FEATS,  to  determine 
temperature  and  thermal  stress  distributions  near  contact  patches  on  the  seal 
ring  surface  [9,10].  It  was  found  tha<-  very  high  temperatures  occur  at  the 
contact  patches  owing  to  frictional  heating.  The  regions  affected  by  these 


high  temperatures  are  quite  small,  but  thermal  stresses  in  those  regions  are 
large.  The  stresseii  are  compressive,  however,  so  cannot  directly  cause  thermo¬ 
cracks  of  the  type  observed  in  our  studies.  It  was  concluded  that  these  large 
compressive  thermal  stresses  cause  substantial  plastic  deformation  near  contact 
patches.  Upon  movement  of  the  contact  patch,  residual  tensile  stress  would 
result  which  could  easily  cause  cracks  like  those  observed,  especially  if  the 
seal  face  material  is  brittle  or  contains  a  brittle  phase  [9]. 

The  analytical  study  resulted  in  a  proposed  thermocrack  mechanism  which 
could  account  for  all  observed  characteristics  of  thermocracks.  The  analysis 
relied,  however,  on  a  crucial  assumption  -  that  the  formation  of  localized 
contact  patches,  or  thermoelastic  instabilities,  precedes  thermocracking.  Al¬ 
though  such  small,  highly  loaded,  contact  patches  had  been  observed  in  various 
laboratory  configurations  by  other  researchers,  there  had  been  no  direct  ob¬ 
servation  of  thermoelastic  instabilities  in  face  seals.  The  only  evidence  had 
been  ex-post- facto,  as  in  the  early  phase  of  this  work. 
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EXPERIMENTAL  STUDY  OF  CONTACT  PATCHES  IN  FmCE  SEALS 

Although  the  hot  spots  resulting  from  thermoelastic  instability  have  been 
studied  analytically  by  a  number  of  researchers  [6],  there  has  been  little  direct 
evidence  of  their  existence  in  seals.  The  most  dramatic  experimental  verification 
of  thermoelastic  instabilities,  by  Dow  and  Stockwell  [12],  was  for  a  case  of  a 
thin  blade  loaded  against  the  outside  of  a  rotating  drum.  The  blade  was  oriented 
parallel  to  the  axis  of  the  drum  so  the  drum  surface  was  moving  perpendicular 
to  the  blade.  The  results  of  that  study  are  not  directly  applicable  to  face 
seals  because  the  kinematics  are  different  from  those  in  a  face  seal.  Banerjee 
and  Burton  [12]  conducted  tests  on  a  ring  on  disk  configuration  kinematically 
similar  to  a  face  seal,  and  observed  the  formation  of  thermoelastic  instabilities. 
The  materials  used  in  their  tests,  aluminum  ring  on  glass  disk,  are  different 
from  those  used  in  face  seals,  however,  so  again  the  direct  applicability  of 
their  results  to  face  seals  is  open  to  question.  Netzel  [8]’  has  reported  evidence 
of  hot  spot  formation  due  to  thermoelastic  instability  in  a  variety  of  face  seals. 
His  observations  were  made  on  seals  which  had  been  disassembled  after  suffering 
from  excessive  wear,  unacceptable  leakage,  and/or  high  frictional  power  loss. 

The  observations  provided  much  evidence  that  instabilities  had  occurred  during 
operation  of  the  seals  and  that  the  resulting- hot  spots  had  contributed  to  the 
seal  failure.  They  did  not  provide,  however,  the  final  proof  of  the  presence 
of  thermoelastic  instabilities  in  operating  face  seals  -  observation  of  the  hot 
spots  while  they  were  occurring. 

The  purpose  of  this  phase  of  the  study  was  the  development  of  an  experi¬ 
mental  technique  to  observe  the  hot  contact  spots  resulting  from  thermoelastic 
instability  during  operation  of  actual  mechanical  face  seals.  The  technique 
could  then  be  used  to  determine  the  size  of  the  spots,  the  motion  of  the  spots 
relative  to  each  of  the  two  contacting  seal  faces,  and  the  influence  of  various 
operating  and  material  parameters  on  hot  spot  characteristics.  It  was  also 
desired  to  devise  and  implement  a  means  for  measuring  temperature  of  the  con¬ 
tacting  surfaces  within  the  hot  spots. 
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EXPERIMENTAL  APPARATUS 

A  monitor  was  desired  which  could  determine  the  number,  size,  location, 
and  relative  motion  of  patches  of  contact  on  the  interface  of  an  actual  face 
seal  during  operation.  It  was  also  desired  that  the  monitor  be  easy  to  use 
and  not  modify  the  operating  characteristics  of  the  face  seal.  Temperature 
measurement  within  the  contact  patch  was  also  desired. 

Although  a  variety  of  contact  probes  are  available,  none  were  deemed 
suitable  for  this  application  and  it  was  instead  decided  to  design  and  construct 
one.  Recall  that  the  major  characteristic  of  thermoelastic  instabilities  is 
nonuni form  contact.  It  follows  that  the  areas  of  contact  and  non-contact,  if 
suitably  different,  would  act  like  the  opening  and  closing  of  an  electrical 
switch  if  current  were  selectively  passed  through  a  narrow  interfacial  region 
between  primary  and  mating  seal  rings,  provided  that  one  of  the  rings  is  con¬ 
ductive.  Implementation  of  this  idea  is  described  schematically  in  Figure  1. 

A  fine  wire  is  implanted  in,  but  insulated  from,  the  stationary  seal  ring 
(carbon  in  this  case).  The  probe  is  then  levelled  flush  with  this  seal  surface 
through  polishing  techniques.  A  graphite  electrical  brush  is  placed  against 
the  outer  diameter  of  the  rotating  mating  ring  to  inject  current  into  the  system. 
Current  flows  from  the  brush  through  the  ring  to  the  probe  wire  only  when  there 
is  direct  contact  at  the  probe  location.  If  small  thermal  asperities,  resulting 
from  thermoelastic  instability,  are  present  on  the  mating  ring  surface,  current 
will  flow  through  the  system  only  when  the  asperities  are  in  contact  with  the 
probe  wire.  An  intermittent  output  (current  flow)  will  thus  be  produced 
(Figure  1).  If  the  contact  patches  are  stationary  with  respect  to  the  moving 
seal  ring,  the  output  will  be  periodic  in  nature,  repeating  once  each  revolution. 
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The  seals  used  in  this  test  program  were  commercial ly-avail able  mechani¬ 
cal  face  seals  designed  for  turbine  engine  applications.  The  seal  has  a  carbon 
graphite  primary  ring  and  a  mating  ring  made  of  440C  stainless  steel.  The 
sealing  surface  has  a  diameter  of  approximately  5  cm.  and  a  width  of  about 
2.5  mm.  To  implant  the  contact  probe,  a  .22  mm  hole  was  drilled  through  the 
carbon  ring  and  a  .178  mm  diameter  constantan  wire,  surrounded  by  a  thin 
coating  of  epoxy  insulation,  was  embedded  in  the  hole  with  epoxy.  The  sealing 
surface  of  the  ring  was  ground,  lapped  and  polished  to  a  surface  finish  of  less 
than  .3ym  to  insure  that  the  end  of  the  wire  was  flush  with  the  surface  and 
that  the  surface  was  at  least  as  flat  as  a  new  commercial  seal  face. 

It  should  be  noted  that  although  the  above  description  was  for  insertion 
of  the  contact  probe  in  the  carbon  ring  to  monitor  contact  between  it  and 
contact  patches  on  the  moving  metallic  mating  ring,  a  similar  procedure  can  be, 
and  indeed  was,  followed  in  which  the  rings  were  switched.  That  is,  a  contact 
probe  was  embedded  in  the  stationary  metallic  ring  to  monitor  contact  patches 
on  the  moving  carbon  ring.  In  this  way  a  complete  mapping  of  contact  patch 
sizes  and  their  relative  motion  relative  to  each  of  the  seal  rings  could  be 
obtained.  This  was  possible  in  thi s  case  owing  to  the  good  electrical  conduc¬ 
tivity  of  both  ring  materials. 

An  electrical  brush  was  designed  to  supply  the  rotating  seal  ring  with  a 
1  milliamp  current  at  5  volts.  The  graphite  brush  and  its  support  rod  were 
spring  loaded  to  provide  a  constant  contact  force  between  the  brush  and  the 
outer  diameter  of  the  rotating  seal  ring.  Both  seal  rings  were  supported  in 
specially  designed  holders,  each  of  which  could  accommodate  either  carbon  or 
metallic  ring.  Each  holder  contained  a  plexiglass  insulating  disk  to  prevent 
current  leakage  from  metallic  components  of  the  holders  to  ground.  The  lower 
holder  also  contains  a  means  for  introduction  of  pressurized  fluid  to  the  in- 


terior  of  the  seal.  Figure  2  shows  a  carbon  primary  seal  ring,  with  contact 
probe,  in  position  in  the  lower  holder*. 

The  upper  holder,  containing  the  rotating  seal  ring,  was  attached  to  the 
spindle  of  a  drill  press.  The  drill  press  allowed  variation  of  the  rotational 
speed  of  the  moving  ring  in  8  speeds  ranging  from  15.7  /sec  to  188.5  /sec, 
(surface  velocity  of  moving  ring  ranged  from  .393  "*/sec  to  4.71  "’/sec).  The 
drill  press  also  allowed  control  of  the  contact  pressure  at  the  seal  interface 
by  means  of  axial  force  application.  Although  the  spindle  speed  remained 
approximately  constant  during  a  test,  a  small  variation  in  speed  would  produce 
a  drift  in  the  output  signal  which  could  cause  errors  in  the  determination  of 
contact  patch  location.  To  overcome  this,  a  simple  tachometer  was  attached  to 
the  spindle  to  trigger  the  oscilloscope  trace.  The  tachometer  consisted  of  a 
notched  circular  plate  attached  to  the  spindle  shaft  and  rotating  between  a 
stationary  LED  light  source  and  a  light  receptor.  Once  each  revolution  the 
rotating  notch  allows  the  passage  of  photons  to  the  receptor  and  this  generates 
a  pulse  to  trigger  the  oscilloscope.  Since  the  trigger  always  occurs  at  the 
same  position  on  the  ring  circumference,  the  output  signal  remains  stationary 
with  respect  to  the  rotating  seal  ring.  Any  movement  of  the  output  pulses 
therefore  represents  a  movement  of  the  contact  patches  with  respect  to  the 
moving  ring  and  that  movement  can  be  continuously  monitored. 

The  problem  of  measuring  the  seal  surface  temperature  within  the  contact 
patches  was  solved  by  using  a  dynamic  thermocouple.  A  thermocouple  is  formed 
by  the  junction  of  two  dissimilar  metals,  with  the  voltage  difference  between 
the  metals  being  linearly  related  to  the  temperature  difference  between  the 
hot  junction  and  the  other,  cold,  junction  of  the  thermocouple  circuit.  In  a 
dynamic  thermocouple  one  of  the  two  metals  at  the  hot  junction  is  moving  rela¬ 
tive  to  the  other.  In  this  case  the  constantan  probe  wire  served  as  one  side 
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of  the  thermocouple  circuit  while  the  metallic  mating  ring,  provided  it  was 
an  iron-or-copper-based  material,  couTd  serve  as  the  other.  With  the  voltage 
source  through  the  graphite  brush  turned  off,  contact  between  a  thermal 
asperity  on  the  moving  metallic  ring  and  the  probe  wire  constitutes  a  thermo¬ 
couple  junction  and  the  small  resulting  voltage  between  mating  ring  and  con- 
stantan  wire  is  proportional  to  the  junction  temperature.  The  junction  voltage 
was  measured  using  a  potentiometer-amplifier  circuit.  Shielding  of  all  wires 
was  used  to  limit  signal  noise,  but  the  thermocouple  output  signal,  being  of 
low  input  voltage,  was  still  of  lower  quality  than  that  from  the  contact  probe. 
The  thermocouple  was  calibrated  by  monitoring  its  static  response  to  several 
known  temperatures  (0®  and  100®C). 

A  test  chamber  was  built  around  the  lower  seal  holder  to  collect  leakage 
of  the  sealed  fluid.  At  the  base  of  the  chamber  was  a  ball  and  seat  which 
maintained  alignment  between  primary  and  mating  rings.  The  entire  lower  part 
of  the  assembly  -  stationary  seal  ring  and  holder,  test  chamber,  and  ball  and 
seat  system  -  rests  upon  a  thrust  bearing  to  allow  free  rotational  motion.  The 
rotation  is  impeded  by  two  strain-gaged  cantilever  arms  fixed  to  the  drill  press 
platform.  The  output  from  the  strain  gages  was  calibrated  to  enable  determina¬ 
tion  of  frictional  torque.  The  nonnal  force  applied  to  the  seal  faces  was 
applied  by  hanging  known  weights  from  a  pulley  attached  to  the- drill  spindle 
lowering  mechanism.  From  this  known  normal  force  and  the  measured  frictional 
torque,  a  friction  coefficient  could  be  determined. 

A  schematic  diagram  of  the  entire  test  apparatus  is  shown  in  Figure  3. 

More  details  of  the  test  system  can  be  found  in  Reference  14. 

TEST  PROCEDURES  AND  MATERIALS 

Three  different  types  of  tests  were  run  in  this  program: 
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Figure  3  Schematic  diagram  of  test  apparatus 
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(1)  Rotation  of  three  different  metallic  mating  rings  against  a  station¬ 
ary  carbon  primary  ring  at  three  different  test  speeds  in  air  under  dry,  un¬ 
pressurized  conditions. 

(2)  Rotation  of  the  carbon  ring  against  a  stationary  mating  ring  under 
dry,  unpressurized  conditions. 

(3)  Rotation  of  a  metallic  ring  against  the  stationary  carbon  ring  with 
pressurized  water  as  a  sealed  fluid. 

The  first  series  of  tests  used  mating  rings  made  fr'om  440C  stainless 
steel  (the  production  material),  hardened  beryllium  copper,  and  52100  bearing 
steel.  The  properties  of  these  materials  are  given  in  Table  1.  Each  of  the 
rings  was  rotated  at  three  different  speeds:  62.8  '^®‘*/sec,  125.7  '^^‘^/sec,  and 
188.5  ’^®‘^/sec. 

To  answer  the  question  of  the  ring  surface  upon  which  the  thermal  asperi¬ 
ties  resided,  the  carbon  ring  was  also  rotated.  These  tests  were  done  with  a 
stationary  metallic  ring  {440C  stainless  steel)  and  with  the  carbon  ring  rotating 
at  four  different  speeds,  the  three  listed  above  plus  31.4  /sec.  In  this  case 
the  contact  probe  was  embedded  in  the  stainless  steel  ring,  but  Insulated  from  it. 

The  only  tests  with  pressurized  water  were  run  with  a  rotating  mating  ring 
m^ue  from  440C  stainless  steel  and  a  stationary  carbon  ring.  The  rotational 
speed  was  125.7  /sec.  Conventional  tap  water  at  a  pressure  of  .35  MPa  was 
Introduced  at  the  interior  of  the  face  seal  after  the  seal  rings  were  in  contact 
but  before  rotation  began. 

The  test  procedure  was  as  follows.  Prior  to  a  test  both  mating  and  primary 
seal  ring  were  lapped,  polished  to  a  surface  roughness  of  .25Mm  or  better,  and 
cleaned  with  acetone.  Proper  functioning  of  the  contact  probe  was  insured  by 
microscopic  examination  and  an  electrical  check,  and  the  temperature  circuitry 
was  checked  for  electrical  bias.  Torque  measuring  circuitry  was  also  checked. 


The  drill  spindle  was  lowered  and  the  desired  normal  force  (65N  in  most  cases) 
was  applied.  The  drill  press  was  turned  on,  with  the  rotational  speed  having 
previously  been  set  to  the  desired  value.  With  the  tachometer/ trigger  system 
on,  monitoring  of  the  contact  probe  output  on  the  oscilloscope  began.  At 
desired  time  intervals,  photographs  of  the  output  trace  were  taken  and  the 
then  current  source  was  turned  off  so  that  the  thermocouple  output  could  be 
monitored.  Photographs  of  that  output  signal  were  also  taken.  The  tests  con¬ 
tinued  until  either  the  signals  remained  unchanged  for  at  least  15  minutes  or 
until  the  contact  probe  stopped  functioning  (see  below).  Most  tests  lasted  at 
least  30  minutes. 

RESULTS  AND  DISCUSSION  OF  EXPERIMENTAL  STUDY 

Rotating  mating  ring  -  dry.  The  first  tests,  with  a  stainless  steel  mating 
ring  rotating  at  188.5  ""^^/sec,  were  aimed  at  proving  the  viability  of  the  con¬ 
tact  probe,  which  was  located  in  the  stationary  carbon  ring.  As  noted  earlier, 
the  mating  ring  had  been  lapped  and  polished  prior  to  testing.  Despite  this, 
the  contact  probe  displayed  evidence  of  distinct  contact  patches  iitiuediately 
upon  operation  (with  no  sealed  fluid).  In  the  first  test  there  was  inmediate 
indication  of  two  distinct  contact  spots,  one  well-defined  and  the  other,  about 
120*  around  the  circumference,  smaller  in  size.  Soon  thereafter,  a  third  spot, 
about  120®  from  each  of  the  first  two,  became  visible.  These  three  spots  grew 
slightly  in  size  until  they  reached,  after  about  15  minutes  of  testing,  the 
size  indicated  in  Figure  4(a).  The  figure  shows  that  the  output  from  the 
contact  probe  is  approximately  a  squarewave  form  with  three  peaks,  100-120® 
apart,  which  repeat  with  each  revolution  of  the  mating  ring.  This  indicates 
that  three  patches  of  contact  are  present  on,  and  rotating  with,  the  mating 
ring  surface. 
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Figure  4a  Figure  4b 

Comparison  between  contact  probe  output  (4a,  left)  and  appearance  of  mating 
ring  surface  (4b,  right). 

440C  Stainless  steel  mating  ring,  188.5  rad/sec,  dry. 
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Comparison  between  contact  probe  output  (Fig.  5a,  left)  and  mating  ring  surface 
appearance  (Fig.  5b,  right). 

440C  Stainless  steel  mating  ring,  188.5  rad/sec,  dry. 
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The  test  was  stopped,  and  the  seal  disassembled,  just  after  the  photo 
of  Figure  4(a)  was  taken.  Observation  of  the  mating  ring  surface  showed  that 
there  were  three  distinct  wear  patches  on  the  surface  (see  Figure  4tb)).  Note 
the  very  close  correlation  between  the  size  and  distribution  indicated  by  the 
contact  probe  and  that  visible  on  the  ring  surface. 

A  second  test,  under  the  same  conditions,  but  with  freshly  polished  seal 
surfaces,  yielded  the  results  shown  in  Figure  5.  Again  there  were  two  small 
spots  apparent  immediately  after  rotation  began,  but  this  time  they  were  180® 
apart  and  those  two  spots  grew  to  the  steady  state  size  shown  in  Figure  5(a). 
The  mating  ring  surface,  upon  disassembly  of  the  seal  soon  after  Figure  5(a) 
was  taken,  showed  evidence  of  contact  in  patches  of  about  the  same  size  and 
in  the  same  location  as  those  indicated  by  the  contact  probe.  This  test,  and 
that  shown  in  Figure  4,  show  conclusively  that  contact  conditions  in  a  face 
seal,  as  determined  either  by  the  contact  probe  or  by  seal  disassembly  and 
observation,  are  not  uniform,  but  that  the  number  of  patches  of  actual  contact 
cannot  be  uniquely  predetermined. 

Although  the  temperature  probe  (dynamic  thermocouple)  was  not  in  operation 
in  the  tests  described  above,  it  was  used  successfully  in  numerous  other  tests. 
Results  of  one  such  test  are  shown  in  Figure  6.  It  can  be  seen  that  the  regions 
of  temperature  indication  coincide  with  the  contact  patch  locations  Cno  thermo¬ 
couple  junction  is  formed  where  contact  does  not  exist),  and  that  the  tempera¬ 
ture  generally  reaches  a  maximum  near  the  trailing  edge  of  the  contact  patch. 
This  is  in  agreement  with  theoretical  predictions  of  surface  temperature  in 
fast  moving  contacts  in  seals  [9]. 

An  extended  series  of  tests  was  run  with  three  different  mating  ring 
materials,  and  at  three  different  mating  ring  speeds.  In  all  tests,  no  matter 
what  the  mating  ring  material  or  the  sliding  speed,  the  contact  probe  showed 


Figure  6  Comparison  between  contact  probe  output 
(above),  and  thermocouple  output  (below) 
Test  E,  440C  stainless  steel  mating  ring 
125.7  rad/sec,  dry. 


19 


evidence  of  small,  distinct  contact  spots  being  present  as  soon  as  seal  oper¬ 
ation  began.  This  occurred  in  spite  of  all  our  efforts  to  lap  and  polish  the 
seals  before  operation  and  it  is  probably  due  to  a  very  slight  initial  waviness 
of  the  seal  faces.  Very  soon  after  the  test  began,  however,  the  spot  size 
began  to  grow,  both  by  expansion  of  the  patches  and  by  formation  and  joining 
together  of  new  contact  spots.  A  steady  state  patch  size,  or  at  least  one  in 
which  no  significant  growth  occurred  during  a  15  minute  period,  was  reached 
within  30-60  minutes  for  all  seals.  In  each  case  these  steady  state  patches 
appeared  to  be  rotating  with  the  mating  ring,  showing  no  motion  with  respect 
to  that  ring.  The  growth  of  asperities  from  their  initial  size  (within  1  min¬ 
ute)  to  their  steady  state  size  (after  about  30  minutes)  is  depicted  in  Figure  7. 

A  summary  of  the  characteristics  of  the  steady  state  contact  patches  for 
each  of  the  test  cases  is  given  in  Table  2,  with  typical  results  for  each  of 
the  seal  materials  being  given  in  Figures  7(b),  8,  and  9.  A  study  of  thft 
results  shows  that  although  there  is  no  unique  set  of  patch  characteristics 
associated  with  any  set  of  test  parameters,  there  are  some  significant  trends. 

For  a  given  seal  velocity,  the  52100  bearing  steel  ring  tends  to  have  the 
smallest  patch  lengths,  the  highest  patch  temperatures  and  generally  the  most 
patches  of  contact,  while  the  beryllium  copper  ring  tends  to  have  the  largest 
and  coolest  patches.  Complete  contact  over  the  entire  circumference  (157  mm) 
was  obtained  in  a  test  at  125.7  "^^^/sec  with  the  beryllium  copper  ring,  but 
only  at  the  lower  speed  (62.8  '"®‘^/sec)  with  the  440C  stainless  steel  ring,  and 
never  with  the  ring  made  from  52100  bearing  steel.  In  general,  higher  patch 
temperatures  accompany  higher  velocities  and/or  smaller  spot  sizes.  The  friction 
coefficient  was  found  to  be  approximately  the  same  for  all  materials  and  veloc¬ 
ities  in  this  series  of  tests,  ranging  from  .10  to  .12.  The  total  heat  generation 
rate  is,  therefo; e,  approximately  proportional  to  velocity.  The  measured  patch 


Characteristics  of  Steady- State  Contact  Patches  for  Various  Mating  Ring  Materials 
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Figure  7  Sequence  of  output  traces  from  contact  probe, 
showing  growth  of  contact  patches  to  steady 
state  size.  Surface  temperatures,  shown  above 
contact  probe  output,  decrease  as  patch  sizes 
Increase. 

Stainless  steel  mating  ring,  188.5  rad/sec. 
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Figure  9  Contact  probe  output  (steady  state) 

52100  bearing  steel  mating  ring,  188.5  rad/sec,  test  0. 


temperatures  are  not  directly  porportional  to  the  velocity  though,  because 
at  a  given  total  heat  generation  rate  temperature  is  increased  by  a  decrease 
in  spot  size  or  by  a  decrease  in  velocity  [9]. 

In  general,  it  can  be  stated  that  the  long-term  presence  of  thermal 
asperities  (or  thermoelastic  instabilities)  is  suppressed  more  by  beryllium 
copper  than  by  stainless  steel  and  least  with  52100  alloy  steel.  The  major 
influences  promoting  the  growth  of  thermal  asperities  are  heat  conduction  and 
wear.  In  Table  1  it  can  be  seen  that  beryllium  copper  is  the  most  conductive 
and  softest  of  the  three  materials.  The  higher  conductivity  of  beryllium 
copper  aids  in  the  dissipation  of  frictional  heat,  leading  to  lower  hot  spot 
temperatures  than  with  the  other  materials.  This  in  turn  leads  to  less  thermal 
dilatation  and  less  driving  force  for  thermal  asperity  formation.  The  lower 
hardness  of  beryllium  copper  leads  to  lower  wear  resistance  and,  therefore, 
easier  removal  of  thermal  asperities  once  they  have  formed,  52100  bearing  steel, 
the  most  wear  resistant  of  the  three  materials,  tends  to  have  thermal  asperities 
which  are  resistant  to  growth  by  wear.  That  material  also  has  the  highest 
product  of  modulus  of  elasticity  and  thermal  expansion  coefficient,  a  factor 
which  determines  the  stress  and  deformation  response  of  a  material  to  a  tempera- 
tu»*£  rise. 

Burton  [6]  has  stated  that  patch  length  is  proportional  to  thermal  con¬ 
ductivity  and  inversely  proportional  to  sliding  velocity,  friction  coefficient, 
and  the  product  of  modulus  of  elasticity  and  thermal  expansion  coefficient,  and 
that  patch  length  is  decreased  by  increased  wear  resistance.  Although  the 
parametric  influences  in  our  tests  are  similar  to  those  predicted  by  Burton, 
attempts  to  use  his  equation  [6],  which  did  not  include  a  wear  coefficient 
term,  to  predict  our  measured  contact  lengths  proved  fruitless.  The  sizes 
of  thermal  asperities  predicted  by  the  Burton  equation  were  much  smaller  than 
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those  measured  experimentally.  In  addition,  that  equation  predicted  that 
52100  bearing  steel  has  a  larger  patch  size  at  a  given  velocity  than  44QC 
stainless  steel,  in  direct  opposition  to  our  experimental  observations. 

Perhaps  these  inaccuracies  could  be  remedied  by  the  inclusion  of  a  determinable 
wear  resistance  term  in  Burton's  equation. 

Rotating  carbon  primary  ring  -  dry.  No  motion  of  the  contact  patches  relative 
to  the  metallic  mating  ring  was  noted  in  the  first  test  series,  in  which  the 
mating  ring  was  rotating.  In  this  set  of  tests  a  contact  probe  was  placed 
In  the  stationary  mating  ring  and  the  carbon  primary  ring  was  rotated.  The 
aim  was  to  see  if,  under  those  different  kinematic  conditions,  there  was  any 
motion  of  the  contact  patches  relative  to  the  metallic  mating  ring.  The  re¬ 
sults  were  negative.  Either  an  indication  of  contact  was  observed  at  the 
contact  probe,  in  which  case  it  remained  in  contact  throughout  the  30  minute 
test,  or  no  contact  was  observed  at  the  probe  location  throughout  the  test. 
Similar  results  were  found  at  all  velocities,  which  ranged  from  31.4  to  188.5 
'^®^/sec.  Upon  disassembly  of  the  seal  after  the  tests,  however,  there  were 
worn  areas  on  the  mating  ring  surfaces,  similar  to  those  shown  in  Figures  4(b) 
and  5(b),  in  many  cases.  Thus,  there  was  evidence  that  thermal  asperities  had 
been  present  on  the  mating  ring  surface,  sometimes  at  the  probe  location  and 
sometimes  elsewhere,  but  the  bumps  did  not  seem  to  be  moving  on  the  mating  ring. 

The  fact  that  these  tests  produced  no  good  evidence  that  thermal  asperities 
were  present  on,  and  rotating  with,  the  carbon  ring  surface,  is  probably  due  to 
the  relatively  low  wear  resistance  of  the  carbon  material.  Any  such  thermal 
bumps  which  may  have  formed  on  that  surface  were  probably  worn  away  nearly 
immediately. 

Rotating  metallic  ring  with  sealed  liquid.  Most  face  seals  are  used  to  seal  a 
liquid,  not  dry  air  as  in  the  earlier  tests,  so  the  evidence  of  contact  patch 


formation  presented  earlier  does  not  necessarily  apply  to  those  liquid  lubri¬ 
cated  seals.  For  that  reason  a  series  of  tests  was  run  using  the  same  test 
apparatus  but  with  pressurized  water  being  present,  as  a  sealed  fluid,  at  the 
inside  diameter  of  the  seal.  A  contact  probe  was  implanted  in  a  stationary 
carbon  primary  ring,  while  a  440C  stainless  steel  mating  ring  was  rotated. 

The  seal  faces  were  prepared,  and  the  tests  run,  in  the  same  manner  as  in  the 
dry  tests. 

Results  of  these  liquid  lubricated  tests  were  much  less  satisfactory  than 
those  of  the  dry  tests.  During  seal  operations,  water  was  present  as  a  very 
thin  lubricating  film  between  the  seal  faces.  This  water,  which  was  slightly 
ionized  and  thereby  conductive,  enabled  current  to  flow  between  the  contact 
probe  and  the  mating  ring  even  when  no  contact  occurred  at  the  probe  location. 
This  current  leakage  produced  a  continuous  output  signal  of  nearly  5  volts, 
and  this  just  about  obliterated  output  signals  from  the  contact  probe.  There 
was,  however,  some  evidence  that  signals  from  contact  patches  were  superimposed 
on  the  leakage  signal.  The  signals  suggest  that  the  fluid  film  is  continually 
being  broken  by  small  patches  of  solid-to-solid  contact.  The  contact  patches 
appeared  to  be  transitory  in  nature,  however,  in  contrast  to  the  long  lasting, 
stable  patches  which  were  observed  in  dry  contact.  It  is  possible  that  the 
liquid  film  aids  in  the  cooling  of  the  hot  contact  patches,  leading  to  a  rela¬ 
tively  short-lived  existence  of  the  patches. 

Work  is  now  in  progress  to  overcome  the  current  leakage  problems  en¬ 
countered  with  liquid  lubricated  seals  and  to  expand  these  preliminary  results. 
It  should  be  mentioned  that  somewhat  similar,  but  unrelated,  current  leakage 
problems  occurred  in  several  dry  tests,  forcing  the  early  curtailment  of  the 
test.  In  those  cases  the  leakage  was  caused  by  wear  debris  from  the  carbon 
seal  ring  which  formed  a  bridge  from  the  constantan  probe  wire  across  the 


surrounding  insulation  to  the  carbon  seal  ring  in  which  the  wire  was  embedded. 
In  those  few  cases  the  problem  was  remedied  by  a  repolishing  of  the  seal  sur¬ 
faces. 

CONCLUSIONS  FROM  EXPERIMENTAL  STUDY 

1.  A  contact  probe/thermocouple  apparatus  has  been  built  which  enables  the 
monitoring  of  contact  conditions  on  the  sealing  surfaces  of  mechanical  face 
seals  during  seal  operation.  The  probe  system  was  applied  successfully  to 
measure  the  size,  distribution,  and  surface  temperature  of  patches  of  solid-to- 
solid  contact  on  the  seal  faces. 

2.  Contact  patches  resulting  from  thermoelastic  instabilities  have  been  ob¬ 
served  to  exist  on  the  surface  of  operating  mechanical  face  seals  under  both 
dry  and  liquid  lubricated  conditions.  At  least  in  dry  operation,  the  patches 
tend  to  remain  stationary  with  respect  to  the  metallic  mating  ring  of  the 
seal,  whether  that  ring  be  stationary  or  rotating. 

3.  Contact  patches  observed  on  beryllium  copper  mating  rings  tend  to  be 
larger,  fewer  in  number,  and  at  a  lower  temperature  than  those  found  with  mating 
rings  made  from  440C  stainless  steel  or  52100  bearing  steel.  The  52100  rings 
iiad  the  smallest  and  hottest  contact  patches  at  all  seal  velocities. 

4.  Based  on  the  results  of  these  tests  it  appears  that  the  mating  ring 
materials  most  conducive  to  the  formation  of  small,  hot,  long  lasting  contact 
patches  (or  thermoelastic  instabilities)  are  those  with  low  thermal  conductivity 
high  wear  resistance,  and  high  modulus  of  elasticity  and  coefficient  of  thermal 
expansion.  A  combination  of  heat  conduction  and  wear  can  lead  to  the  enlarge¬ 
ment  of  the  contact  patches,  resulting  in  more  uniform  contact  and  lower  contact 
temperatures.  A  decrease  in  velocity  also  results  in  more  uniform  contact, 
larger  contact  patches,  and  lower  contact  temperatures. 


ANALYTICAL  STUDY  OF  TEMPERATURES  NEAR  CONTACT  PATCHES 

The  goal  of  this  analytical  study  was  the  determination  of  temperature 
distributions  in  mechanical  face  seals.  It  was  earlier  [9,10]  shown  that  these 
temperatures  cause  thermal  stresses  which  are  the  probable  cause  of  thermo¬ 
cracking  of  the  seal  faces.  It  was  desired  to  determine  the  influence  of  various 
operating  and  material  parameters  on  seal  temperatures.  To  accomplish  these 
goals,  a  thermal  analysis  of  the  contacting  rings  was  carried  out  using  finite 
element  methods. 

Based  on  the  results  of  the  experimental  study,  it  was  assumed  that  thermo¬ 
elastic  instability  precedes  thermomechanical  failure  (e.g.,  thermocracking), 
and  that  contact  is  concentrated  at  three  (3)  patches  located  at  120“  from  each 
other  on  the  contact  interface.  It  was  also  assumed  that  the  patches  are 
stationary  with  respect  to  the  rotating  mating  ring,  as  was  observed  in  our 
experimental  program.  The  case  of  contact  patches  remaining  stationary  with 
respect  to  the  carbon  primary  ring  had  been  studied  earlier  [9,10].  The  patches 
were  estimated  from  microscopic  observations  of  seal  surfaces  [9]  to  be  about  .1 
to  1  mm  long  and  dimensions  in  that  range  were  used  in  this  study. 

It  was  assumed  that  each  of  the  three  contact  patches  encountered  identical 
contact  conditions,  so  only  one  120“  sector  of  the  seal  needed  to  be  analyzed 
(see  Figure  10).  A  finite  element  model  of  the  shaft  seal  is  shown  in  Figure  11. 
Overall  dimensions  and  thermal  boundary  conditions  were  modelled  after  those  of 
a  typical  marine  shaft  seal.  The  finite  element  mesh  was  two-dimensional  (axial 
and  circumferential  directions),  with  variable  thickness  of  the  elements  being 
used  to  account  for  heat  flow  in  the  radial  direction.  Elements  as  small  as 
10  pm  xp30  m  were  used  in  the  contact  zone,  with  elements  as  large  as  24  mm 
X  15  mm  used  elsewhere  in  the  mesh.  The  mesh  used  in  the  model  contained  694 
nodes  and  644  elements. 

For  this  model  a  uniformly  distributed  frictional  heat  source  was  assumed 
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to  be  located  at  the  contact  interface.  Convection  boundary  conditions  were 


assumed  to  occur  along  the  non-contacting  portions  of  the  seal  interface,  and 
fixed  temperatures  were  applied  to  the  boundary  segments  representing  the  back 
sides  of  the  seal  ring  inserts.  Material  properties  represented  those  of  actual 
seal  ring  materials,  and  those  properties  were  varied  for  the  purpose  of  deter¬ 
mining  their  influence  on  surface  temperatures.  Typical  properties  and  boundary 
conditions  are  shown  in  Figure  11.  A  full  parametric  study,  which  involved 
variation  of  most  of  those  properties,  was  performed  to  provide  guidance  for 
future  seal  material  selection. 


FINITE  ELEMENT  THERMAL  ANALYSIS  PROGRAM  (THERMAP) 

The  finite  element  program  used  in  this  study  was  based  on  techniques 
developed  earlier  for  studying  surface  temperatures  in  sliding  systems  [11]. 

The  program  has  recently  been  modified  and  expanded,  and  its  capabilities  have 
been  proven  by  comparison  with  closed  form  solutions  for  a  variety  of  problems  [15], 
This  new  program,  THERMAP,  enables  the  transient,  steady  state,  or  quasi-steady 
state  analysis  of  temperature  distributions  in  moving  or  stationary  bodies,  or 
in  systems  composed  of  both  moving  and  stationary  components.  In  a  transient 
analysis,  boundary  conditions,  including  convection,  surface  heat  flux,  or  pre¬ 
scribed  temperature  may  vary  with  time,  as  may  velocity  and  internal  heat  genera¬ 
tion  rate.  Material  properties  may  be  a  function  of  temperature. 

The  equations  on  which  THERMAP  is  based  are  developed  in  Appendix  A.  A 
flow  chart  for  the  program  is  presented  in  Appendix  B,  and  user  instructions 
for  the  program  are  in  Appendix  C. 

The  program's  performance  was  verified  by  application  to  ten  different 
problems  for  which  closed- form  solutions  were  available  [15].  Each  problem  was 
designed  to  prove  out  a  different  capability  of  THERMAP.  Excellent  agreement 
was  attained  between  theory  and  numerical  prediction  in  all  cases  [15].  It  was 
noted  that  numerical  inaccuracies  may  occur  in  the  solution  of  problems  involving 
either  high  moving  body  velocities  or  excessively  large  time  increments  (for 
transient  cases).  The  inaccuracies  for  both  situations  were  characterized  and 
criteria  for  their  avoidance  were  developed  [15].  Since  the  first  of  these 
two  inaccuracies  may  occur  in  seal  analyses,  that  problem  was  studied  in  some 
detail  and  a  discussion  of  the  problem  and  its  solution  follows  later  in  this 
report. 

One  additional  capability  which  was  added  to  THERMAP  so  that  systems  like 

those  in  Figures  10  and  11  could  be  studied  was  "node  tying."  This  procedure 
enables  the  equating  of  unspecified  temperatures  at  any  two  (or  more)  nodes  in 

the  finite  element  mesh.  For  example,  the  seal  shown  in  Figure  10  has  three 


B INDICATES  ASSUMED  CONTACT 


»  PATCH  LOCATIONS 

- DESIGNATES  PORTION  OF  SEAL 

COMPONENTS  MODELLED  IN 
FINITE  ELEMENT  STUDY 


Figure  10,  Sketch  of  primary  seal  ring,  showing  section  studied  in  thermal  and 
stress  analyses. 
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contact  patches  on  its  surface.  Each  patch  is  considered  identical,  and  the 
end  of  one  120®  sector  has  the  same  temperature  as  the  beginning  of  the  next 
sector.  The  temperature  distribution  is,  therefore,  periodic  in  space.  In 
order  to  analyze  only  one  120®  sector,  as  in  Figure  11,  the  nodes  of  one  end 
of  the  seal  sector  must  be  "tied"  to  those  at  the  other  end  so  that  the  tempera¬ 
tures  at  the  two  ends  will  be  equal. 


FACTORS  INFLUENCING  SURFACE  TEMPERATURES  IN  FACE  SEALS 

Using  the  THERMAP  program  described  above,  temperatures  were  calculated 
for  the  face  seal  model  shown  in  Figure  11.  A  typical  temperature  distribution, 
calculated  for  a  contact  width  2b=.06  cm,  had  isotherms  shown  in  Figure  12. 

In  the  figure  the  mating  ring,  which  moves  from  left  to  right,  is  at  the  top, 
while  the  stationary  primary  ring  is  at  the  bottom.  The  contact  patch  is 
moving  with  the  mating  ring.  It  can  be  seen  that  the  maximum  surface  tempera¬ 
ture  occurs  near  the  rear  (trailing  edge)  of  the  contact.  Temperature  gradients 
in  the  depth  direction  are  greater  in  the  primary  ring,  with  respect  to  which 
the  contact  patch  is  moving.  Much  of  the  frictional  heat  is  therefore  being 
conducted  into  the  primary  ring.  The  high  temperatures  and  high  gradients  in 
both  seal  rings  lead  to  large  thermal  stresses  which  are  predominantly  compressive 
in  nature.  In  the  earlier  phase  of  this  project  it  was  concluded  that  the  large 
compressive  stresses  created  near  the  hot  patches  induce  localized  plastic 
yielding  on  the  mating  ring  seal  face  [9,10].  As  the  hot  patches  move  to  another 
spot  on  the  seal  -face,  owing  to  wear,  the  thermal  stress  is  released  and 
residual  stresses  occur  in  the  plastically  deformed  region.  It  was  hypothesized 
that  these  tensile  residual  stresses  are  responsible  for  thermocrack  initiation 
[9,10]. 

To  investigate  the  influence  of  the  modelling  assumptions  and  to  observe 
the  relative  importance  of  the  various  material  properties,  a  parametric  study 
of  all  system  variables  was  carried  out.  The  results  of  39  different  calculations 
are  shown  in  Table  3.  From  runs  1,  21,  and  22  or  5  and  6,  it  can  be  seen  that 
surface  temperatures  are  approximately  proportional  to  total  frictional  heat 
generation  rate  (=  friction  force  x  velocity).  At  a  given  heat  generation  rate, 
temperatures  decrease  with  increasing  velocity  (runs  1-5),  but  the  influence 
of  velocity  is  not  as  strong  as  that  of  heat  flux.  Thus,  at  a  constant  friction 
force  the  temperature  increases  as  the  velocity  increases  (runs  1  and  6  or  3  and 
22),  This  is  in  agreement  with  our  experimental  results  (Table  2).  One  conclusion 
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of  this  study  is  that  a  reduction  in  seal  velocity  would  reduce  seal  temperatures 
and,  consequently,  thermal  stresses.  This  could  be  accomplished  by  a  decrease 
of  seal  diameter. 

The  study  showed  that  at  a  given  heat  generation  rate  and  velocity,  a 
reduction  in  contact  patch  size  would  result  in  higher  patch  temperatures 
(runs  23-27).  This  also  agrees  with  our  experimental  results  (Figure  9,  for 
example).  The  influence  of  contact  width  on  temperature  is  quite  a  bit  smaller 
than  that  of  total  heat  generation  rate,  so  heat  flux  (heat  generation/unit  area) 
is  less  of  a  factor  than  total  heat  generation.  Convective  cooling  was  found 
to  have  a  relatively  small  influence  on  surface  temperature  (runs  1,  18-20) 

4 

with  no  significant  cooling  occurring  until  convection  coefficient  h  becomes 
very  high.  Ambient  temperature  Ta  has  little  influence  over  a  reasonable  range 
of  values  (runs  1,  15-17).  Assumed  prescribed  temperature  boundary  conditions 
have  a  bit  larger  influence  than  convection  so  that  bulk  cooling  could  help 
lower  surface  temperatures.  The  change  would  not  be  dramatic,  however  (runs  1, 
7-14). 

Material  properties,  especially  those  of  the  primary  seal  ring,  have  a 
substantial  influence  on  surface  temperatures.  An  increase  in  thermal  conduc¬ 
tivity  of  the  primary  ring  would  cause  a  reduction  in  surface  temperatures 
(runs  1,  33-35),  while  an  increase  in  that  material's  thermal  diffusivity, 
which  normally  accompanies  an  increase  in  conductivity,  would  have  a  lesser, 
but  opposite,  effect  (runs  1,  36-39).  The  diffusivity  of  the  mating  ring 
material,  on  which  the  hot  spot  is  stationary,  has  no  effect  on  surface  temper¬ 
atures  (runs  31-32),  but  the  conductivity  of  that  material  is  a  significant 
factor  (runs  1,  28-30).  This  dependence  of  surface  temperatures  on  mating  ring 
conductivity  was  also  noted  in  our  experimental  study  (Tables  1,  2).  Consider¬ 
ing  these  results,  along  with  the'  interrelationship  between  thermal  conductivity 
and  thermal  diffusivity,  the  best  approach  to  reducing  surface  temperatures  in 
seals  might  be  to  choose  a  higher  conductivity  mating  ring  material.  An  im- 
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provement  in  the  primary  ring  conductivity  would  not  be  quite  as  effective 
unless  the  conductivity  increase  could  be  accomplished  without  a  corresponding 
increase  in  diffusivity. 


ABLE  3*  ParafflatrLc  Study  for  KechanicaX  Face  Seal  Analysis 

(Subscript  p  refers  to  properties  of  primary  seal  ring  material 
and  subscript  m  refers  to  pr«^erties  of  mating  ring  material 


EFFECT  OF  MOVING  SURFACE  VELOCITY  ON  THERMAL  ANALYSIS 
Evidence  of  Numerical  Instabilities  at  High  Velocity 

The  equations  solved  in  the  finite  element  analysis  of  surface  temperature 
are  derived  from  Fourier's  law  for  heat  conduction,  equation  (1)  in  Appendix  A, 

At  low  surface  velocities,  the  contribution  of  the  term  VVT  in  that  equation  is 
quite  small,  the  governing  differential  is  approximately  parabolic,  and  the 
resulting  finite  element  equations  (equations  (12)  in  Appendix  A)  are  approxi¬ 
mately  symmetric.  Few  problems  result  in  the  solution  of  the  equations  and  the 
results  are  quite  accurate  as  long  as  the  mesh  is  fine  enough  to  model  tempera¬ 
ture  gradients  and  the  time  increment  is  small  enough  to  follow  transient  changes 
[15].  When  the  velocity  is  large,  however,  the  first  derivative  terms  become 
significant,  the  governing  equation  is  to  some  extent  hyperbolic,  and  the  re¬ 
sulting  finite  element  equations  are  nonsyiranetric.  It  has  been  found  by  numerous 
investigators  in  recent  years  [e.g.,  18],  that  numerical  oscillations  may  develop 
in  the  solution  of  equations  similar  to  these  at  high  velocities.  Such  equations 
often  called  convective-diffusion  equations  because  of  the  presence  of  the 
convective  or  advective  term  V*VT,  appear  in  a  variety  of  transport  problems, 
and  numerical  oscillations  are  encountered  in  both  finite  difference  and  finite 
element  solutions  of  the  equations  in  fluid  flow  and  heat  transfer  investigations 
In  order  to  study  the  magnitude  of  the  problem,  a  simple  sliding  contact 
problem  (metallic  slider  on  metallic  flat)  was  analyzed  under  quasi-steady  state 
conditions  using  the  finite  element  model  shown  in  Figure  13a  on  the  left.  A 

blowup  of  the  mesh  at  the  contact  zone  is  shown  in  Figure  13a  on  the  right. 

Heat  generation  due  to  friction  was  modeled  by  a  uniform  heat  flux  within  the 
contact  zone.  A  fixed  temperature  of  22*C  was  set  at  the  lower  left  corner  of 

the  mesh.  Figures  20b  -  20c  show  the  isotherms  predicted  by  THERMAP  for  values 

of  velocity  ranging  from  V=0  to  V=100,000  cm/s.  Isotherms  are  plotted  both  in 
the  .complete  mesh  ancl’in  the  contact  region  at  regular  intervals  of  AT.  The 
Isotherms  drawn  tn  Figures  20b  -  20e  seem  quite  reasonable  in  both  sliding  and 
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stationary  bodies.  For  these  low  velocities  the  isotherms  and  temperature 
magnitudes  are  affected  by  velocity,  but  no  solution  difficulties  were  en¬ 
countered.  When  the  velocity  reaches  50  cm/s  for  this  problem,  however,  spatial 
oscillations  begin  to  appear  in  the  moving  body  temperature  field  and  at  slightly 
higher  velocities  (Figs.  20g-20k)  the  surface  temperature  predictions  are 
definitely  inaccurate.  At  even  higher  velocities  (Figure  20il-20m)  the  tempera¬ 
ture  magnitudes  seem  reasonable,  but  spatial  oscillations  reappear  in  the  moving 
body.  Thus  it  became  apparent  that  the  application  of  THERMAP  to  problems  of 
high  sliding  velocity  could  lead  to  invalid  temperature  predictions. 

To  study  the  problem  in  more  detail,  another  sample  problem  was  posed.  In 
this  problem,  the  model  of  which  is  shown  in  Figure  14a,  a  uniformly  distributed 
heat  flux  q  was  applied  to  a  band  on  the  surface  of  a  moving,  two-dimensional 
slab.  Temperatures  in  the  body  were  analyzed  for  a  variety  of  different  conditions 
and  the  results  of  some  of  the  runs  are  given  in  Table  4. 

A  study  of  the  data  in  Table  4  shows  that  oscillations  may  occur  whenever 
the  velocity  of  the  moving  body  relative  to  the  heat  source  reaches  a  critical 
value.  That  critical  value  is  a  function  of  the  diffusivity  of  the  material  of 
that  body  and  of  the  size  of  the  finite  elements.  In  particular,  it  appears 
that  oscillations  occur  whenever  the  element  Peclet  number  is  equal  to  or  greater 
than  2,  Pg  *  ^  1  ^ 

The  critical  Peclet  number  appears  to  be  applicable  to  the  largest  element 
in  a  region  where  there  are  significant  temperature  gradients.  In  runs  18-20, 
for  example,  it  was  found  that  oscillations  first  began,  not  at  the  heat  source, 
but  a  slight  distance  away  where  the  elements  were  larger  (and  P^  was  larger). 

If  the  velocity  continues  to  increase,  oscillations  could  move  into  regions 
where  the  elements  are  smaller  when  the  critical  Peclet  number  for  those  elements 
is  reached.  This  critical  value  of  2  for  the  element  Peclet  number  is  the  same 
number  determined  by  other  researchers  [19]  for  analysis  of  convective  diffusion 
equations  using  linear  elements.  That  study  also  showed  that  the  critical 
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Velocity  effect  in  frictionally  heat 
generating  contact  problem.  (System  shown  on  left,  contact  zone  on  rt) 
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value  may  be  slightly  different  for  different  element  interpolation  functions 
[19].  The  results  of  our  study  (e.g.,  run  11)  also  show  that  the  critical 
Peclet  number  value  may  be  even  lower  than  2  if  the  temperature  gradients  within 
the  elements  are  large,  owing  to  high  heat  flux,  for  example. 

It  can  also  be  seen  that  results  at  Peclet  numbers  above  the  critical  value 
may  not  be  too  inaccurate  unless  the  oscillations  are  occurring  in  the  region 
of  interest,  usually  the  heat  flux  zone.  As  was  the  case  in  the  earlier  study 
(Figure  13)  at  Peclet  numbers  much  higher  than  2,  Run  17  for  example,  the  results 
can  be  quite  reasonable  and  relatively  oscillation-free.  The  worst  results  seem 
to  occur  at  Peclet  numbers  between  4  and  40.  There  are  a  great  many  potential 
applications  in  that  region.  In  such  cases,  the  problem  of  numerical  oscillations 
and  inaccuracies  could  be  a  serious  deterrent  to  the  use  of  finite  element  thermal 
analysis  programs  like  THERMAP  in  surface  temperature  analysis. 

There  are  two  possible  remedies  to  this  problem.  The  simplest  is  to  make 
the  finite  element  mesh  smaller,  at  least  in  the  region  where  numerical  oscilla-  ' 

I 

tions  are  occurring.  This,  in  fact,  may  be  the  best  solution  in  many  cases  be-  ! 
cause  it  may  enable  more  accurate  representation  of  temperature  gradients  in  the 
region  of  interest.  A  large  reduction  in  mesh  size  may,  however,  lead  to  an  | 

excessive  increase  in  computer  time  and  memory  requirements.  An  alternative  is 
to  modify  the  numerical  procedure  in  order  to  avoid  solution  oscillations.  This 
approach  has  been  the  focus  of  intensive  study  recently. 

Techniques  for  Handling  Convective-Diffusion  Equations 

The  problem  of  numerical  oscillations  at  high  Peclet  number  occurs  in  both 
finite  element  and  finite  difference  solution  of  equations  of  the  convective- 
diffusion  type.  Such  equations  occur  in  a  variety  of  transport  problems,  in¬ 
cluding  fluid  flow,  seepage,  and  heat  conduction.  Because  of  this,  numerous 
Investigators  have  been  attempting  in  recent  years  to  find  ways  to  overcome  the 
oscillation  problems.  In  this  work  we  are  most  concerned  with  techniques  for 


TABLE  4 

Results  of  Thermal  Analysis  of  Moving  Slab  Subject  to  Distributed 

Heat  Source 


Heat 

Flux 

Run  q  AXi 

AXp 

Diffu- 

sivity 

D 

Veloc¬ 

ity 

V 

Ml  T 

D  max 

1  100  .1 

.1 

1 

10 

1  111.5 

2  100  .1 

.1 

1 

20 

2  '108.2 

3  100  .1 

.1 

1 

40 

4  106.0 

4  100  .1 

.1 

.667 

10 

1.5  114.1 

5.  100  .1 

.1 

.5 

10 

2  116.4 

0 

0 

.1 

.25 

10 

4  ;124.1 

( 

7’  100  ^U5 

.05 

1 

35 

1.75  104.4 

! 

8'  100  i.05  1  .05 

1 

J 

40 

2  il04.1 

9,  100  |.05 

.05 

1 

50 

2.5  103.7 

1 

10; 1000  i.05 

.05 

: 

'  ! 

30 

1.5  147.0 

‘ 

11:1000  ' .05 

4  , 

.05 

1  , 

35 

1.75  143.5 

12:1000  ;.05 

.05 

’  j 

40 

2  140.9 

t  j 

13,1000  .05 

.05 

1 

100 

5  jl27.1 

14:1000  .05- 

.05 

1 

200 

10  120.3 

1 

15  1000  ,.05 

.05 

1 

400 

20  119.7 

16  1000  i.05 

.05 

1 

800 

40  120.9 

n'lOOO  |.05 

.05 

1 

1600 

80  119.2 

18  1000  .1 

.05 

1 

20 

2  156.5 

19  1000  j.l 

.05 

1 

25 

2.5  150.9 

20; 1000  j.l 

.05 

1 

30 

3  146.6 

. 

21;  1000  !.l 

.05 

1 

40 

4  140.3 

_ Comments _ 

good  results 

slight  subsurface  oscillations 
oscillations  on  and  beneath  surface 
good  results 

slight  subsurface  oscillations 
oscillations  on  and  beneath  surface 
good  results 

beginnings  of  subsurface  oscillations 
subsurface  oscillations 
good  results 

beginnings  of  subsurface  oscillations 
subsurface  oscillations 
oscillations  under  contact 
slight  oscillations  under  contact 
no  subsurface  oscil.,  slight  surf.rippli 
no  subsurface  oscil.,  ripple  on  surface 

II  II  II  It  II  II 

oscillation  just  beginning(in  AXi  regipfj 
slight  oscil.  in  AX]  region 
oscil.  beginning  to  touch  AX2  region 
oscillations  under  contact  (AX2  region) 
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Figure  14 


Finite  Element  Test  Models  for  Study  of  Velocity  Effect 


(a)  2-Dimensional  Model  (Moving  Slab) 

Heat  Flux  q  moving  at  Velocity  V  relative  to  slab 


V 


(b)  l-Dimens1onal  Model  (Moving  Rod) 

Heat  Flux  q  moving  at  Velocity  V  relative  to  Rod 
Same  mesh  dimensions  as  In  (a) 


use  with  finite  element  methods  and,  for  the  moment,  in  two-dimensional  studies 
with  bilinear  finite  elements. 

A  variety  of  techniques  have  been  developed  for  use  in  this  type  of  situ¬ 
ation  in  the  past  few  years.  Among  them  are  the  following: 

1.  The  method  of  Finn  and  Varoglu  [20],  which  employs  space-time  finite 
elements  and  incorporates  the  characteristics  of  the  hyperbolic  portion  of  the 
partial  differential  equation. 

2.  A  method  developed  by  O'Neill  and  Lynch  [21],  in  which  the  heat  con¬ 
duction  equations  are  formulated  in  terms  of  a  moving  coordinate  system.  Although 
the  convection  or  advection  term  is  absent  in  such  a  formulation,  the  finite 
element  mesh  must  be  continuously  deforming. 

3.  A  procedure  recently  proposed  by  Kanarochos  [22],  in  which  the  con¬ 
ventional  Galerkin  finite  element  equations  (eqns.  12  in  Appendix  A),  are  solved, 
but  special  boundary- layer  type  elements  are  used  in  that  portion  of  the  moving 
body  in  which  the  temperature  gradients  are  significant. 

4.  Upwinding  finite  element  schemes,  used  by  Neinrich  et.al.  [23],  and 
others.  In  that  technique  the  finite  element  equations  are  developed  using  a 
weighted  residual  method  with  a  weighting  function.  Wj  different  from  the 
interpolation  functions  Ni.  The  weighting  functions  used  by  Heinrich,  et.  al. 
[23],  were  of  the  form  Wi  =  N-f  +  a  F^*  where  a  is  a  function  of  the  element 
Peclet  number. 

5.  The  streamline  upwind  formulation  developed  by  Brooks  and  Hughes  [25], 
in  which  the  weighting  functions  (see  #4  above)  are  constants  which  act  only  in 
the  velocity  direction. 

6.  A  different  type  of  upwinding  developed  by  Baliga  and  Patankar  [24], 
in  which  both  weighting  function  and  interpolation  function  vary  with  Peclet 
number.  In  their  formulation  the  functions  vary  exponentially  in  the  velocity 
direction,  but  become  equal  to  the  linear  interpolation  functions  Ni  when  the 
velocity  decreases  to  zero. 


A  study  of  each  of  these  techniques  and  their  applicability  to  heat  con^ 
duction  problems  with  moving  heat  sources  revealed  that  the  upwinding  techniques 


are  the  easiest  of  the  above  procedures  to  implement  and  possibly  the  most 
suitable  for  this  application.  '  Two  of  the  techniques  (#5  and  6)  were  implemented 
in  THERMAP. 

Application  of  upwinding  technique  to  THERMAP 

In  the  original  formulation  of  the  finite  element  equations  in  THERMAP  (see 
Appendix  A),  the  Galerkin  procedure  was  employed  with  weighting  functions  Ni, 
which  are  also  the  interpolation  functions  used  in  the  approximation  of  the 
temperature  across  an  element.  For  the  linear  triangle  element,  the  interpola¬ 
tion  functions  are  given  by  [16]: 

Nf  =  ^  [(XJ  VXk  Xj)  +  X(Yj-V|,)  +  Y(X|<-yj)]  (2) 

where  (X,*,  Y^),  (Xj,  Yj),  and  (X^,  Y^)  are  the  coordinates  of  the  of 

the  element's  three  corner  nodes 
and  A  is  the  area  of  the  element. 

In  our  version  of  the  Brooks-Hughes  technique,  the  interpolation  functions 

remain  the  same  but  the  weighting  functions  become: 

3N.  3N.  5 

Wi  =  Ni  +  IT  (Vx  jy!-  +  Vy  -jl  )/V‘=  (3) 

where  Vx  and  Vy  are  the  x  and  y  components  of  velocity 


and  IT  is  a  constant  which  depends  on  the  element  Peclet  number  [25] 
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.025 
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2 
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near  source 

2 
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1 

5 
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.1 

1 

2.5 

98.89 

100.82 
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.05 

1 

1 

.1 

105.47 

105.47 

Baliga  good  results 

.02 

1 

5 
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100.82 
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m 

.01 

1 

.1 

99.73 

100.07 

"  slight  osc.  near 

source 

.05 

1 

1 

.1 

105.47 

105.47 

Brooks  good  results 

.01 

1 

10 

100.07 

100.07 

II  H  H 

.05 

1 

1000 

100.00 

100.00 

II  II  11 

TABLE  6 


Moving  Slab  with.  Distributed  Heat  Source 
-  Results  of  Analysis  with  Upwinding 


Velocity 

V 

V  Xi 

Tempera  1 

ture  at  Node  66  (~  ^max 

D 

No  Upwinding 

Baliga 

Brooks 

0 

0 

605.39 

605.26 

1  605.39 

1 

.001 

.0001 

604.91 

604.94 

'  604.91 

.01 

.001 

602.64 

602.61 

602.64 

.1 

.01 

580.73 

580.58 

580.68 

1 

.1 

421.51 

421.48 

.  421.38 

10 

1 

180.05 

179.91 

i  179.84 

20 

2 

156.5  * 

156.77 

156.53 

40 

4 

140.3  ** 

140.5  * 

139.85 

100 

10 

123.21  ** 

125.34  ** 

123.37 

1000 

100 

1 

- 

104.37 

♦  Oscillations  just  about  beginning. 
**  Oscillations  in  solution. 


When  these  weighting  functions  are  used  In  a  weighted  residual  procedure 
(the  Petron-Galerkin  procedure),  finite  element  equations  are  obtained  which 
are  Identical  to  those  given  In  Appendix  A  except  for  a  modified  Kv^.  term 
(In  equation  3,  Appendix  A).  The  modified  term  reduces  to  the  Kv^.j  given  In 
Appendix  A  for  the  case  of  zero  velocity.  For  non-zero  velocities,  the  in¬ 
clusion  of  PI  weights  the  resulting  stiffness  matrix  contribution  In  the  direction 
of  positive  velocity. 

In  the  Baliga-Patankar  method,  the  weighting  function  is  set  equal  to  the 
Interpolation  function,  but  both  are  different  from  the  relations  given  in 
eqn.  (2).  A  new  coordinate  system  X,Y  is  defined  in  the  x,y  plane  with  the 
X  axis  oriented  in  the  velocity  direction.  A  new  variable  Z  is  defined  [24] 
such  that: 

Z  =  -§-  exp  (X-Xmax)  -  1]  (4) 

where  D  is  the  material's  diffusivity 

and  Xmax  is  the  maximum  (most  positive)  X  coordinate  location 
of  the  element's  nodes. 

Using  this,  the  interpolation  function  becomes: 

"tj  '  \  -  Zk  »j)  *  Z(Yj  -  Y^)  *  Y(Z|^  -  Z.)]  (5) 

where  Z^. ,  Z^,  Z|^  are  the  Z  values  at  the  element's  nodes 

and  is  the  area  of  the  triangle  formed  by  • 

(Z^,  Y^),  (Z2m  Y2),  (Z3,  Y3) 

This  interpolation  function  is  then  used  in  eqn.  (3)  of  Appendix  A  to  get  the 
finite  element  equations  for  this  formulation.  Because  of  the  non-linear  nature 
of  Z,  the  integrals  for  each  element  must  be  evaluated  numerically. 

Modified  versions  of  THERMAP  were  written  to  include  the  above  modifications 
in  the  finite  element  equations.  The  programs,  along  with  the  original  version, 
were  applied  to  the  two  test  problems  shown  in  Figure  14.  One  of  the  problems 
is  the  same  2-dimensional  problem  studied  earlier  (Table  4),  while  the  other 


is  a  problem  involving  1-dimensional  heat  flow  in  a  moving  rod  (Figure  14b). 

The  results  of  some  of  the  tests  are  presented  in  Tables  5  and  6. 

In  Table  5  the  predictions  of  the  original  version  of  THERMAP  are  compared 
with  those  of  the  Brooks-Hughes  aiid  Baliga-Patankar  versions.  Also  included 
are  analytical  predictions  from  a  closed  form  solution  available  for  this  case 
[26].  All  three  versions  of  the  program  gave  good  predictions  for  element 
Peclet  numbers  less  than  2,  in  accordance  with  equation  (1).  When  the  Peclet 
number  of  the  largest  element  in  the  mesh  reached  or  exceeded  2,  however, 
oscillations  began  to  appear  in  the  numerical  results  using  the  original  THERMAP 
version.  The  oscillations  were  accompanied  by  inaccuracies  in  the  temperature 
predictions  for  nodes  just  ahead  of  the  heat  source.  The  oscillations  and 
inaccuracies  occurred  whether  the  Peclet  number  was  increased  by  an  increase 
in  element  size,  an  increase  in  velocity,  or  a  decrease  in  diffusivity.  The 
Baliga-Patankar  version  gave  good  predictions  for  element  Peclet  numbers  less 
than  about  8,  a  four-fold  improvement  over  the  original  version.  At  higher 
Peclet  numbers,  however,  the  results  deteriorated,  producing  the  same  type  of 
Inaccuracies  and  oscillations  noted  at  lower  Peclet  numbers  with  the  original 
version  (e.g..  Run  8).  The  Brooks-Hughes  version,  on  the  other  hand,  gave  good 
predictions  over  the  entire  range  of  Peclet  numbers  from  0  to  10,000.  No 
oscillations  were  noted  in  any  of  the  solutions  with  that  technique  and  all 
results  agreed  well  with  theoretical  predictions. 

Results  for  the  two  dimensional  problem  (Fig.  14a)  are  given  in  Table  6. 

Results  are  somewhat  less  satisfying  than  for  the  one-dimensional  heat  flow 

problem.  As  was  seen  earlier,  the  original  version  of  THERMAP  produced  good 

results  for  element  Peclet  numbers  less  than  2.  With  the  Baliga-Patankar  version, 

oscillation-free  results  continued  until  Peclet  number  reached  4  but  oscillations 

occurred  in  solutions  at  higher  Peclet  numbers.  Again  the  Brooks-Hughes  method 

gave  the  best  results,  with  oscillations  not  beginning  until  the  largest  element 
* 

Peclet  number  reached  about  60.  Some  oscillations  occurred  at  higher  Peclet 


numbers,  but  in  all  cases  they  were  less  severe  than  those  found  with  the 
other  techniques  at  the  same  Peclet  number.  Slight  differences  were  noted 
between  the  predictions  of  the  different  programs,  but  the  differences  were 
not  great.  No  closed- form  solution  exists  with  which  to  compare  the  predictions. 


CONCLUSIONS  FROM  ANALYTICAL  STUDY 

The  finite  element  surface  temperature  analysis  program  THERMAP,  developed 
and  used  in  this  study,  has  proven  to  give  good  predictions  of  seal  temperature 
distributions  for  a  wide  variety  of  conditions.  Use  of  the  program  in  a 
parametric  study  showed  that  seal  surface  temperature  is  approximately  propor¬ 
tioned  to  total  heat  generation  rate  (=  fixation  force  X  seal  surface  velocity). 

A  reduction  in  seal  velocity,  by  a  reduction  in  seal  diameter,  would  result  in 
a  temperature  decrease.  Smaller  decreases  in  temperature  could  be  accomplished 
by  an  increase  in  convective  cooling  or  by  bulk  cooling  of  the  seal  rings. 
Material  properties  of  the  seal  rings  were  found  to  have  a  significant  influence 
on  seal  surface  temperatures,  with  higher  thermal  conductivity  of  both  primary 
and  mating  seal  rings  leading  to  substantial  decreases  in  temperatures.  If.  the 
increase  in  primary  seal  ring  conductivity  is  accompanied  by  a  corresponding 
increase  in  that  material's  diffusivity,  however,  much  of  the  decrease  in  surface 
temperature  would  be  wiped  out. 

The  thermal  analysis  showed  that  numerical  oscillations  and  inaccuracies 
can  occur  if  the  element  Peclet  number  of  the  primary  seal  model  exceeds  2  in 
the  vicinity  of  significant  temperature  gradients.  This  can  occur  with  either 
high  surface  velocities,  low  thermal  diffusivity,  or  large  element  sizes.  The 
problem  is  similar  to  that  encountered  by  other  researchers  in  recent  years  with 
equations  of  the  convective-diffusion  type.  Several  methods  were  studies  as  a 
means  to  avoid  the  numerical  oscillations  and  it  was  found  that  the  streamline 
upward  procedure  developed  by  Brooks  and  Hughes  leads  to  a  substantial  improve¬ 
ment  in  the  solution  at  high  Peclet  numbers  without  excessive  computing  effort. 
That  procedure  has  now  been  incorporated  in  THERMAP. 
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APPENDIX  A 


THEORETICAL  DEVELOEflENT  OF  EQUATIONS 

Consider  Fourier's  lav  for  heat  conduction  in  an  anisotropic 
solid  moving  with  velocity  V: 

V  •  KVT  ♦  Q  -  /»C(^T/ht)  -  V  •  VT  «  0  (1 ) 

Kennedy[tl]  developed  the  finite  element  equations  for  the  quasi¬ 
stationary  condition,  C^T/it)  >  0,  reducing  equation  (l)  to: 

V  KVT  +  Q  -  V  •  VT  -  0  (2) 

The  temperature  solution  must  satisfy  the  following  boundary 
conditions: 

T  *■  Ts  prescribed  temperature  in  specified  regions,  SI 
qs*n  «  KvT«n  prescribed  heat  flux  normal  to  surface,  S2 
b(Ta-T)«n  «  KvT*n  convection  to  Ta  across  surface,  S2 

The  continuum  being  modeled  is  discretized  into  finite  elements 
with  a  temperature,  Ti,  being  associated  at  each  vertex  (node). 
Assuming  that  the  temperature  field  is  linearly  distributed 
throughout  each  element  temperature  is  defined: 

T  -  ^HiTi®  -  [n]®It)® 

<•1 

where  m  is  the  total  number  of  nodes  in  an  element  and  Ni  are  the 
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interpolation  functions[l6].  Applying  the  Galerkin  niethod[t6]  and 
integrating  the  resulting  equation  by  parts,  Kennedy[lO  arrives 
at  the  finite  element  equations  in  matrix  notation: 

([Kt]*[Kh]-[Kv])  {t}  -  {Q}-*-{ql  +  {KTa}  (5) 

where: 

■  1  <i‘[Kx(‘aii/»x)(SNj/di)  +  Ky(9Ni/»y)(»Nj/dy)]dxdy 


Kh  -  I  hNiNjdS 

'v 

..  -/ocj  dNi[VxONj/9x)  ♦  Vy(h»j/iy)]dxdy 


Kv 


Q.  ■  1  dQNidxdy 


■  I  qsNidS 


"•i-  f 


hTaNidS 


for  2<-dimensional  heat  flow  in  the  xy  plane. 

To  solve  3-dimenslnal  axisymmetric  systems  (i.e.  2-diBen- 
sional  heat  flow  in  the  rz  plane),  the  matrix  notation  is 
redefined: 


Kt..  -  2»  lr[Kr(»Ni/9r)(»Nj/^r)  ♦  Kz(9Ni/»z)(®Nj/^z)]drdz 


hNlNJdS 


KTy  -  zvcjr! 
"l  - 

'is*' 


NiEVrC-SNj/^r)  ♦  V20Mj/Ba!)]dr<i2 


rQNldrdz 


qsNidS 


f 


KTa,  -  I  MaNidS 


Hotd  that  the  Integral  over  the  element  surface,  Se,  vill  take  on 
different  values  depending  on  whether  the  problem  is  2-dimension¬ 
al  or  3-D  axlsymmetric . 

for  simplicity  equation  (3)  may  be  rewritten: 


[s]{t1  -  (r1  (4) 

where  [s]  ■  [Kt]+[Kh]+[Kv]  and 
Ul  -  {QMq}  +  {KTal. 

The  desired  temperature  distribution  is  arrived  at  by  so-lving  for 
.{t}  in  equation  (4). 


To  modify  THERNAP  for  handling  transient  thermal  analyses, 
the  finite  element  equations  must  be  derived  from  equation  (I) 
instead  of  equation  (2).  This  modification  changes  equation  (4) 
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[s]1t1  ♦  [c]{Ti  -  {Rj  (5) 


f 

where  [c]  -/oC  \dNiNjdxdy  for  2-D  flow  in  xy  plane  or 

[c]  ■  anyiC^rNiNjdi'dz  for  2-D  axisymmetric  flow  in  rz 


plane < 


Equation  (5)  differs  from  equation  (4)  by  the  addition  of  the 
time  derivative  term, 

To  facilitate  solving  equation  (5)  by  computer,  the 
following  finite  difference  time  stepping  scheme[l7]  is  used: 

'v 

iTlt  -  iT)^^  ♦  ♦  (l-e>lTl^J  (6) 

?or  6*«5,  equation  (6)  reduces  to  the  Crank-Nichol8on[l7] 
approximation . 

Bearrange  equation  (6)  to  solve  for  }t}^: 

-  (l/eAt)({T}^-{T^)  -  ((l-e)/0){T}t^t  (T) 

Consider  equation  (5)  at  time,  t,  and  at  an  earlier  time,  t-Bt 
to  produce  the  following  equations: 

[sl^lTii  ♦  C0](T|,  -  1*)^  (8) 


1 
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# 


t 


where  the  matrix  [c]  is  time  invariant. 

Substituting  equation  (7)  into  equation  (8)  and  rearranging  the 
terms  yields  the  equation: 

(Cs]^(i/eAt)[c]){T}^-  ((i-e)/0)[c]{T^  (i/eAt)[c]{T^  {h}^  (io) 

Solving  for  in  equation  (9)  and  substituting  into  equation 

(10)  results  in  the  final  equation: 

([sVi/eAt)[c]){T}  -  ((i/eAt)[c]-((i-e)/e)[s]){T}  +  ((i-e)/e){R}  ♦  {r}. 

^  *  t-ii  tst  t 

(11) 

Transient  temperature  distributions  are  predicted  by  specifying  6 
and  At  as  well  as  the  normal  boundary  conditions  represented  in 
equation  ( 1 1 ) .  The  transient  solution  is  found  by  an  time 
stepping  process,  choosing  At  small  enough  to  attain  numerical 
accuracy  (see  section  V).  To  solve  for  steady  state  temperatures 
by  equation  (11)  would  require  too  many  -time  steps.  Alterna¬ 
tively,  steady  state  temperatures  are  calculated  by  equation  (4). 
Equation  (ll)  may  now  be  rewritten  in  the  form: 

[k]{t}^  -  If}  (12) 

The  total  conductivity  matrix,  [k],  is  an  NxN  non-symmetric 
banded  matrix,  ft  is  the  total  number  of  nodes  in  the  finite 
element  mesh.  The  asymmetry  of  [k]  is  due  to  the  directional 
property  of  the  velocity  vector  which  contributes  to  the  term 
,  [Kv].  Thus  [k]  is  symmetric  only  when  V  -  0.  (p)  is  a  single 
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column  matrix  with  N  total  values . forming  the  right  hand  side  of 
equation  (ll).  A  gaussian  elimination  technique  is  used  to  solve 
equation  (12)  for  the  temperature  vector  {t}^(N  values)  at  time 
tf  incrementing  with  steps  of^t. 
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APPENDIX  B 

THERMAP  FLOWCHART 


VJ5TART  J 

rtLL  ArvriAY' 
lOFEN  ALL  J.-.TA  f  IL^ 
ICALL  I  Ilf  Ci  1' 


999  CALL  i'R 


CALL  LE^TIB 


[INITIALISE  ALL  AHHAY^ 


READ  ALL  INPUT  DATA*  set 


FORIfi  STIFFNESS  XATSIX 

for  A  C3  elements! 
CALL  TSM 
tor—*  elcmentsi 
CALL  LSM 

[Sil 

IFORM  F3IE  FOR  USE  IN 

ITRAN  S*  0— L-ITR  AN  Sa  1 

Jr 


CALCULATE  AND  rRINT  At 


INITIALi; 


FOKX  load  ^»A’aKIa 


r OR  c  1 0  • 


MODIFY* 


gjj 

Hr ’I 


EOLVE  FCR  TE;.:P  DISIkIBUTION 
from  eqm  iKJfTj  =  ir\  ■ 


PRINT  TEXP  DISTRIBUTION  {T 


NtLA0»0-»-NPLAG»l 

I  I 

CO  TO  999  I 
GO  TO  9974—^ 


set  NFLA! 


t 
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APPENDIX  C 


'■J  c 

4  c 
s  c 

^  C - THERMAL  ANALYSIS  PROGRAM - 

7  C 

8  C 

9  C 

10  C  MODIFICATIONS  LAST  MADE  ON  SEPTEMBER  IS. 1982  BY  F.  E.  KENNEDY 

11  C 

12  C 

13  C  THERMAP  SOLVES  2-DIMENSIONAL  AND  3-D  AXISYMMETRIC  HEAT  CONDUCTION 

14  C  PROBLEMS  USING  THE  FINITE  ELEMENT  METHOD. 

15  C 

16  C  THERMAP  ORIGINATED  FROM  SAMPLE  CODE  PRESENTED  BY  K.H.  HUEBNER 

17  C  IN  HIS  textbook:  'THE  FINITE  ELEMENT  METHOD  FOR  ENGINEERS'. 

18  C  THERMAP 'S  CAPABILITIES  HAVE  SEEN  EXPANDED.  MODIFIED  AND  UPGRADED 

19  C  BY  PROFESSOR  FRANCIS  E.  KENNEDY,  JR.  AND  ROGER  P.  GLOVSKY  <M.E.  1982). 

20  C  AT  THAYER  SCHOOL  OF  ENGINEERING,  DARTMOUTH  COLLEGE. 

21  C 

22  C  THERMAP  CAN  SOLVE  TRANSIENT  OR  STEADY  STATE  PROBLEMS  FOR  THE  FOLLOUING 

23  C  CONDITIONS:  FIXED  TEMPERATURE.  CONDUCTION.  HEAT  FLUX  NORMAL  TO  SURFACE 

24  C  HEAT  FLUX  NORMAL  TO  BOUNDARY,  CONVECTION  ACROSS  SURFACE.  CONVECTION 

25  C  ALONG  BOUNDARY.  INTERNAL  HEAT  GENERATION.  AND  NODE  TYING  (EfiUATING 

26  C  TEMPERATURES  AT  ANY  TUO  NODES).  THE  HEAT  CONDUCTION  PROBLEM  MAY  BE 

27  C  2-D  OR  3-D  AXISYMMETRIC  UITH  MOVING  OR  STATIONARY  BODIES. 

28  C  THE  FINITE  ELEMENT  MESH  MAY  CONTAIN  OUADR ILATERAL , TR I ANGULAR . 

29  C  OR  LINEAR  ELEMENTS.  THE  MATERIAL  PROPERTIES  MAY  VARY  FROM  ELEMENT 

30  C  TO  ELEMENT.  AND  MAY  BE  A  FUNCTION  OF  TEMPERATURE. 

31  C 

32  C  THIS  FORTRAN  PROGRAM  REQUIRES  ACCESS  TO  THE  FILE  PARM. 

33  C  ALL  DIMENSIONS  FOR  ARRAYS  ARE  SET  IN  THAT  FILE  AND  ANY  CHANGES 

34  C  IN  THOSE  DIMENSIONS  MUST  BE  MADE  IN  PARM  BEFORE  COMPILING  AND 

35  C  EXECUTING  THERMAP. 

36  C 

37  C  INPUT  TO  THE  PROGRAM  COMES  FROM  14  DATA  FILES.  THE  FILES  AND 

38  C  THEIR  DATA  REQUIREMENTS  ARE  AS  FOLLOWS  (ALL  FILES  ARE  IN  FREE  FORMAT): 

39  C 

40  C  FILE<7)  'ELFILE'  -  CONTAINS  ELEMENT  DATA 

41  C  ONE  LINE  FOR  EACH  ELEMENT: 

42  C  CELEMENT  NO..  N00E1 ,  N00E2.  N0DE3.  N0DE4.  MATERIAL  NO.] 

43  C  NODES  SHOULD  BE  INPUT  COUNTERCLOCKWISE. 

44  C  FOR  TRIANGULAR  CLEMENTS.  SET  NODES  =>  N0DE4. 

45  C  FOR  LINEAR  ELEMENTS.  SET  N0DE2  >  N0DE3  =  N0DE4. 

46  C  FILE<8)  'BCFILE'  -  CONTAINS  PRESCRIBED  TEMPERATURES 

47  C  ONE  LINE  FOR  EACH  PRESCRIBED  NODAL  TEMPERATURE: 

48  C  CB.C. NUMBER.  NODE.  PRESCRIBED  TEMPERATURE] 

49  C  FILE<9)  'IHFILE'  -  CONTAINS  INTERNAL  HEAT  GENERATION  DATA 

50  C  ONE  LINE  FOR  EACH  ELEMENT: 

51  C  CCONOITION  NUMBER.  ELEMENT.  QI] 

52  C  81  >  HEAT/VOLUMC  (HEAT  SOURCE  HAS  POSITIVE  SIGN) 

53  C  FILEdO)  'SHFILE'  -  CONTAINS  PRESCRIBED  HEAT  FLUX  DATA 

54  C  ONE  LINE  FOR  EACH  HEAT  FLUX  BOUNDARY  CONDITION. 

55  C  IF  HEAT  FLUX  IS  NORMAL  TO  BOUNDARY  SEGMENT.  DATA  LINE  ISl 

56  C  CSEGMCNT  NO.  LEFT  NODE.  RIGHT  NODE.  QH] 

57  C  IF  HEAT  FLUX  IS  NORMAL  TO  ELEMENT  SURFACE,  DATA  LINE  ISi 

58  C  CSEGMENT  NO.  ZERO(O).  ELEMENT  NO.  QH] 

59  C  QH  >  HEAT/AREA  (HEAT  SOURCE  HAS  POSITIVE  SIGN) 

60  C  FlLEdl)  'CBFILE'  -  CONTAINS  PRESCRIBED  CONVECTION  DATA 

61  C  ONE  LINE  FOR  EACH  CONVECTION  BOUNDARY  CONDITION. 

62  C  IF  CONVECTION  OCCURS  ALONG  BOUNDARY  SEGMENT.  DATA  LINE  ISl 

63  C  CSEGMENT  NO.  LEFT  NODE.  RIGHT  NODE.  CONVECTION  COEFF.  AMBIENT  TEMP 

44  C  IF  CONVECTION  OCCURS  ACROSS  ELEMENT  SURFACE.  DATA  LINE  ISl 

65  C  CSEGMENT  NO,  ZERO(O).  ELEMENT  NO,  CONVECTION  COEFF,  AMBIENT  TEMP] 


64 


Appendix  C  (continued) 


C  FILE<12)  'aUFILE'  -  CONTAINS  CONTROL  INFORMATION 


C  LINE  1:  LNN.NE.NBC.NHF.NCB.NHG.NNr.NMA.NTHJ 

C  •  LINE  21  ENCASE, NLISr.NCnP.RATIOrNVELl 

C  LINE  3:  CITRANS. rBEQ.OELT. TEND, MINER J 

C  LINE  4:  CNIT,7INIT,THCrAI 

C  LINE  Ss  ENOSC,nBLlNCR,MNFINCR,MViNCR,MCBINERJ 

C  LINE  6:  CITMAX,ERRNAXI 

C  NN  «  NO.  OF  NODES 

C  NE  =  NO.  or  ELEMENTS 

C  NBC  B  NO.  OF  FIXED  lEMPERATURE  UOUNPART  CONDITIONS 

C  NHF  =>  NO.  OF  SEGMENTS  WITH  PRESCRIBED  HEAI  FLUX 

C  NC8  =  NO.  OF  BOUNDART  SEGMENTS  SUBJECT  TO  CONVECTION 

C  NHC  =  NO.  OF  ELEMENTS  UITH  INTERNAL  HEAT  GENERATION 

C  NNT  =  NO.  OF  NODE  TYING  PAIRS 

C  NMA  -  NO.  OF  DIFFERENT  MATERIALS 

C  NTH  =  NO.  OF  DIFFERENT  THICKNESSES 

C  NCASE  3  1  FOR  2-OIMENSIONAL 

C  »  2  FOR  AXISYMMETRIC 

C  NLIST  -  0  IF  ALL  INPUT  DATA  ARE  TO  BE  PRINTED 

C  =  1  IF  ONLY  BOUNDARY  CONDITIONS  ARE  TO  BE  PRINTED 

C  =  2  IF  NO  INPUT  DATA  ARE  TO  BE  PRINTED 

C  NCMP  >  0  FOR  NO  COMPARISON  OF  ELEMENT  ASPECT  RATIOS 

C  =1  1  FOR  PRINT  OUT  OF  LARGE  ELEMENT  ASPECT  RATIOS 

C  RATIO  =  MAXIMUM  DESIRED  CLEMENT  ASPECT  RATIO  (CHOOSE  >=0) 

C  NVEL  =  b  IF  NEITHER  BODY  IS  MOVING 

C  =  -1  IF  ONE  BODY  IS  TRANSLATING  UITH  RESPECT  TO  OTHER 

C  =  1  IF  ONE  BODY  IS  ROTATING  UITH  RESPECT  TO  OTHER 

C  ITRANS  =  0  FOR  STEADY  STATE  SOLUTION 

C  =1  FOR  TRANSIENT  SOLUTION 

C  TBEG  °  TIME  AT  UHICH  TRANSIENT  SOLUTION  UILL  BEGIN 

C  OELT  =  TINE  INCREMENT  FOR  TRANSIENT  SOLUTION 

C  TEND  °  TIME  AT  UHICH  TRANSIENT  SOLUTION  UILL  END 

C  MINCR  =  NUMBER  OF  TIME  INCREMENTS  DESIRED  BETUEEN  PRINTOUTS 

C  NIT  =  1  FOR  READING  INITIAL  TEMPS  FROM  FILE(T5) 

C  »  0  FOR  UNIFORM-  INITIAL  TEMPS 

C  TINIT  =  INITIAL  TEMPERATURE  (UNIFORM) 

C  THETA  °  UEIGHT  FUNCTION  FOR  NEUMARK  RECURSION  RELATION 

C  NOSC  »  1  FOR  CHANGING  BOUNDARY  CONDITIONS 

C  0  FOR  CONSTANT  BOUNDARY  CONDITIONS 

C  MBCINCR>>  NO  OF  TIME  INCREMENTS  BETUEEN  CHANGING  FIXED  B.C. 

C  MHFINCR=  NO  OF  TIME  INCREMENTS  BETWEEN  CHANCING  H.F.  CONDITIONS 

C  MVINCR  =  NO  OF  TIME  INCREMENTS  BETUEEN  CHANCING  VELOCITY  CONDITIONS 

C  MC8INCR=  NO  OF  TIME  INCREMENTS  BETUEEN  CHANGING  CONVECT.  CONDITIONS 

C  ITMAX-MAX  NO  OF  TEMPERATURE  SOLUTIONS  TO  ACHIEVE  PROPERTY 

C  CONVERGENCE 

C  ERRMAX^MAX  ALLOWABLE  ERROR  IN  MATERIAL  PROPERTIES 

C  FILE(13)  ’MAFILE'  -  CONTAINS  MATERIAL  PROPERTY  DATA 
C  ONE  SET  OF  LINES  FOR  EACH  DIFFERENT  MATERIAL 
G  FIRST  LINE  IN  SET:  CMAT, INAT,NUMTC3 

C  NAT»MATERIAL  NUMBER  (REFERENCED  IN  ELEMENT  FILE) 

C  INAT=  1  FOR  MATERIAL  IN  BODY  MOVING  RELATIVE  TO  HEAT  SOURCE 

C  0  FOR  BODY  STATIONARY  RELATIVE  TO  HEAT  SOURCE 

C  NUMTC^NO  OF  TEMPERATURES  AT  UHICH  PROPERTIES  GIVEN  (MAX>8> 

C  REMAINING  LINES  IN  SET  (NUMTC  MORE  LINES) 

C  TEMPER ATURE,X-CONDUCTIVITY,Y-CONDUCTIVITY,DENSITY, SPECIFIC  HEAT 

C  FILEdA)  ’TOFILE'  -  CONTAINS  THICKNESS  DATA 
C  MULTIPLE  LINES  FOR  EACH  DIFFERENT  THICKNESS: 

C  CTHICKNESS,  NO  OF  ELEMENTS  UITH  THAT  THICKNESS3 

C  CLIST  OF  ELEMENTS  WITH  THAT  THICKNESS! 

C  FILEdS)  'TEFILE'  -  CONTAINS  RESULTING  TEMPERATURE  DISTRIBUTION 
C  CTEMP,  TEMP,  TEMP,  TEMP,  TEMP,  TEMP,  TEMP! 

C  FILEdG)  ’RHFILE'  -  CONTAINS  OLD  RIGHT  HAND  SIDE  MATRIX 
C  CRHSO,  RHSO,  RHSO,  RHSO,  RHSO,  RHSO,  RHS03 

C  FILEd7)  'VEFILE'  -  CONTAINS  VELOCITY  OF  MOVING  BODY 
C  COMEGA,XROT, YROT! 

C  OMEGA  ■  ANGULAR  VELOCITY  IF  BODY  IS  ROTATING 

C  -  ZERO  IF  BODY  IS  TRANSLATING 

C  XROT  ■  X  COORDINATE  OF  CENTER  OF  ROTATING  BODY 

C  X  COMPONENT  OF  VELOCITY  IF  TRANSLATING 

C  YROT  ■  Y  COORDINATE  OF  CENTER  OF  ROTATING  BODY 

c  Y  COMPONENT  OF  VELOCITY  IF  TRANSLATING 

C  FlLEdB)  ’NTFILE'  -  CONTAINS  NODE  TYING  DATA 
C  CBC  NO.,  NODE,  NODE! 

C  FILEd?)  'NOFILE'  -  CONTAINS  NODAL  POINT  COORDINATES 
C  ONE  LINE  FOR  EACH  NODE: 
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Appendix  C  (continued) 
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